arXiv:1507.04193v2 [hep-th] 14 Apr 2016 


Integrability and conformal data of the dimer model 


Alexi Morin-Duchesne*, J0rgen Rasmussen^ Philippe Ruelle* 

*Institut de Recherche en Mathematique et Physique 
Universite catholique de Louvain, Louvain-la-Neuve, B-1348, Belgium 

§ School of Mathematics and Physics, University of Queensland 
St Lucia, Brisbane, Queensland 4072, Australia 

alexi .morin-duchesne® uclouvain.be j .rasmussen@uq.edu.au 
philippe. ruelle @ uclouvain. be 


Abstract 

The central charge of the dimer model on the square lattice is still being debated in the literature. 
In this paper, we provide evidence supporting the consistency of a c = —2 description. Using Lieb’s 
transfer matrix and its description in terms of the Temperley-Lieb algebra TL„ at /3 = 0, we provide 
a new solution of the dimer model in terms of the model of critical dense polymers on a tilted 
lattice and offer an understanding of the lattice integrability of the dimer model. The dimer transfer 
matrix is analysed in the scaling limit and the result for Lq — ^ is expressed in terms of fermions. 
Higher Virasoro inodes are likewise constructed as limits of elements of TL n and are found to yield 
a c = —2 realisation of the Virasoro algebra, familiar from fermionic be ghost systems. In this 
realisation, the dimer Fock spaces are shown to decompose, as Virasoro modules, into direct sums of 
Feigin-Fuchs modules, themselves exhibiting reducible yet indecomposable structures. In the scaling 
limit, the eigenvalues of the lattice integrals of motion are found to agree exactly with those of the 
c = — 2 conformal integrals of motion. Consistent with the expression for Lq — ^ obtained from the 
transfer matrix, we also construct higher Virasoro inodes with c = 1 and find that the dimer Fock 
space is completely reducible under their action. However, the transfer matrix is found not to be a 
generating function for the c = 1 integrals of motion. Although this indicates that Lieb’s transfer 
matrix description is incompatible with the c = 1 interpretation, it does not rule out the existence 
of an alternative, c = 1 compatible, transfer matrix description of the dimer model. 

Keywords: dimer model, critical dense polymers, Temperley-Lieb algebra, integrability, conformal 
field theory, conformal integrals of motion, fermionic be ghosts, Feigin-Fuchs modules. 
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1 Introduction 


The classical dimer model, also known as the perfect matching or domino tiling problem, has a rather 
long history. It was originally introduced in 1937 by Fowler and Rushbrook as a model for the adsorption 
of diatomic molecules on a substrate [1]. Its mathematical formulation is extremely simple and asks 
to enumerate the number of ways the edges of a finite graph can be covered by dimers (dominos) so 
that each vertex is covered by exactly one dimer. This is the so-called close-packed dimer model. A 
common generalisation allows for monomers (or vacancies), corresponding to vertices not covered by 
dimers, and leads to the general monomer-dimer problem. Except for a brief mention later in this 
introduction, this paper is exclusively concerned with the close-packed dimer model, formulated on a 
square lattice grid. 

The first important results regarding the statistical properties of the model, such as partition func¬ 
tions and correlation functions, were obtained in the sixties, mostly in the two-dimensional incarnation 
of the model [2-7]. Dating from the same period, the paper [8] by Lieb reformulates and solves the 
dimer model in terms of a transfer matrix; it will play a crucial role in this paper. 

After these pioneering works, the dimer model has been revisited and generalised in a number of 
ways by both the physics and mathematics communities, improving our understanding of the model 
and of its numerous relations with other problems. Recent developments include arctic circle type 
phenomena [9-11], dimers and amoebae [12], correlations of monomers [13-15], the development of 
quantum dimer models [16-18], trimer and more generally p-mer models [19], and double dimer models 
[20,21], A recent review of mathematical aspects of the dimer model can be found in [22]. 

The two-dimensional dimer model belongs to the class of exactly solvable models, since, in most 
instances, it can be solved exactly and explicitly, by various techniques. Often, such models are in fact 
integrable, in a technical sense, where the integrability is related to the presence of algebraic structures 
such as the Yang-Baxter equation, the presence of a spectral parameter in the transfer matrix or, in 
the scaling limit, the existence of an infinite number of conserved charges. It appears that the dimer 
model has not been explicitly demonstrated to be integrable in any of these senses. Our first goal is 
thus to introduce a spectral parameter in the transfer matrix, in the case where the model is defined on 
a strip. The construction relies on two essential ingredients: (i) the fact that the dimer model carries 
a representation of the Temperley-Lieb algebra with fugacity (3 = 0, as shown recently in [23], and 
(ii) the relation, based on Temperley’s correspondence [24], of its transfer matrix T(a) with a family 
of generalised polymer models described by a two-parameter transfer tangle D(u,v ) pertaining to the 
same Temperley-Lieb algebra. 

As indicated, the dimer transfer matrix T(a) depends on a parameter a, but this parameter is not 
spectral: The corresponding family of transfer matrices is not self-commuting. This is contrasted by the 
family of transfer tangles D(u, v ) where u is a spectral parameter (while v is not). The crucial point is 
that we can rewrite the square of T(a ) as the transfer tangle D(u(a), v(a)) in a spin representation for 
particular values of u = u(a) and v = v(a). From this, it directly follows that the one-parameter family 
of matrices D(u,v(a)), labelled by u , commutes with T(a). The coefficients in the series expansion 
in powers of sin (2 u) then yield lattice integrals of motion , a self-connnuting family of operators all 
commuting with T(a), establishing the integrability of the dimer model. Their number is finite for 
finite system size n, grows linearly with n and becomes infinite in the scaling limit. Furthermore, as 
we will show, the connection with the polymer transfer tangles provides a new solution of the dimer 
model. 

This integrability of the dimer model also gives a way to understand the large n expansion of the 
eigenvalues of (the logarithm of) T(a), which here takes the general form [25-27] 

OO 

E = Eig(- ±logT 2 (a)) = nf hulk + / My + ^ ( L1 ) 

p =i 

where /bulk and /bd y are respectively the bulk and boundary free energies. Although the eigenvalues of 
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T(a) are real, they can be positive or negative. It is therefore easier to work with ^ logT 2 (a) instead 
of log T(a), as in (1.1). 

From an exact computation, we will find that the coefficients a p in (1.1) take the form 

2Ptt 2 p~ 1 

a P = (2p)\ ^p- 1 ^ ^p- 1 ’ (-*-2) 

where the (polynomial) functions \ 2 p -1 (a), p ^ 1, are the same for all energy values. The peculiar 
relationship between T 2 (a) and the matrix realisation of the transfer tangle D(u(a),v(a )) shows that 
the hp-i are related to the eigenvalues of lattice integrals of motion. On general grounds [28,29], this 
is exactly what one would expect. Restricting to the conformal part of the spectrum, it is also expected 
that the numbers l^-i become universal in the continuum limit, in the sense that they correspond 
to the eigenvalues of universal (conformal) integrals of motion Iop-i- The first integral of motion is 
simply 1 1 = L 0 — while the higher-order ones are computable polynomials of degree p in the Virasoro 
modes [30]. The 1/n expansion was also obtained in [31] and conjecturally related to the conformal 
integrals of motion. 

In the particular case of the ground-state energy Eq, the first coefficient (1.2) is given by 

<4°’ = vrAi(a) (A 0 - ^) = ~^\i(a)c eS , (1.3) 

where Aq is the conformal weight of the state to which the ground state converges in the scaling limit. 
The factor Ai(a) is usually referred to as the speed of sound, while the combination c e g = c — 24Ao is 
the effective central charge. 

As the previous discussion already assumes implictly, it is widely believed that the dimer model 
on a square lattice 1 is conformally invariant. However, there has been some discussions about how that 
symmetry is realised, and our second main theme concerns this issue. Many works have been devoted 
to various aspects of this question, often with different conclusions. In fact, two principal proposals 
have been put forward: a conformal description with c = —2 and one with c = 1. Let us briefly review 
some of the arguments in favour of one or the other. 

A number of contributions have concentrated on the analysis of finite-size corrections, for different 
boundary conditions and different geometries. This issue is quite subtle because of the fact, known 
from old results by Ferdinand [7], that the dimer model is very sensitive to the parities of the box 
dimensions in which it is defined (as recalled above, it is even more sensitive to the lattice type). 
This peculiarity requires one to be extremely cautious in the analysis of finite-size data, in particular 
when inferring the central charge of the model. The problem is even more acute when the model is 
defined on a rectangle, in which case the free energy contains, at order p = 1, a second universal term 
^logn. Each corner of the rectangle contributes a term of that form that depends not only on the 
central charge [33], but also on the conformal weight of a boundary condition changing field in case the 
boundary conditions on either side of the corner are different [34]. Somewhat surprisingly, these two 
finite-size corrections in the free energy can be consistently interpreted in a conformal scheme based on 
either candidate for the central charge. Namely, those relying on the equivalence with spanning trees 
support c = —2, on a cylinder [35] and on a rectangle [36,37], whereas the use of the height function to 
describe dimer configurations yields c = 1 [38] . As indicated in the case of the rectangular geometry, 
the main difference between the two approaches stems from the presence or absence of appropriate 
changes of boundary conditions, contributing (or not) to the free energy and affecting the way the 
lattice result is to be interpreted. 

The description in terms of spanning trees and its close connection with the Abelian sandpile 
model [39], widely believed to be described by a c = —2 logarithmic conformal field theory [40], 
substantiates the c = —2 interpretation. Combinatorial properties of spanning trees and spanning 

1 In contrast, the dimer model is not conformally invariant on non-bipartite graphs such as the isotropic triangular 

lattice, where correlations decay exponentially [32]. 
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webs have also been considered within the c = —2 framework in [41,42], Another hint is the connection 
between the dimer model and the critical dense polymer model [43] which will be reinforced in this 
paper: As mentioned above, the dimer transfer matrix squared can be written in terms of a transfer 
tangle pertaining to the Temperley-Lieb algebra with fugacity (3 = 0 [23], which is the value canonically 
associated with c = —2. However, a striking result supporting the alternative view is the proof [44] 
that the fluctuations of the height function are described, in the scaling limit, by a massless free scalar 
field, strongly pointing to c = 1. The use of the Kasteleyn matrix and Pfafhan techniques, essential 
tools in the classical solution of the dimer model [2] and still much in use today [15], is closely related 
to free fermionic theories [32, 45] and adds further support to c = 1. The identification of certain 
operators [46], including the insertion of monomers and dimers, as vertex operators in a free boson 
(Gaussian) theory points in the same direction. 

Results from the computation of dimer correlation functions similarly allow for dual represen¬ 
tations. Dimer correlations can be naturally accounted for in a c = —2 conformal theory [40] (the 
operator inserting a dimer can be explicitly written in terms of symplectic fermions), but can also be 
easily interpreted as c = 1 correlators [38]. This is in sharp contrast to the corresponding results for 
monomer correlations. The correlation functions for an arbitrary (even) number of boundary monomers 
have been computed and shown to be equal to the correlators of a complex (Dirac) free fermion [13], 
indicating c = 1. Likewise, the bulk monomer correlators, computed in [14], are consistent with the 
monomer field identification proposed in [46] and do not appear to be naturally accounted for within 
a c = — 2 framework. 

Together, these observations suggest that there exist two conformal descriptions of the dimer 
model: one with c = —2 and one with c = 1. The descriptions would presumably be related to 
different degrees of freedom and come with their own limitations. In this sense, the name dimer model 
can be ambiguous and cover different situations. Depending on whether one keeps track of spanning 
tree connections, includes height degrees of freedom, monomers or still other features, one of the two 
descriptions may impose itself. 2 

Our objective here is to present a self-consistent close-packed dimer model fully compatible with 
a conformal description based on c = —2. Lieb’s transfer matrix is central to our approach. 

Relying on spin degrees of freedom, the transfer matrix can be written as a spin chain, in terms of 
Pauli matrices, and diagonalised by standard techniques. The conformal spectrum can then be easily 
computed, but the spectrum, or even the characters of the representations involved in the model, is 
in general not enough to identify the central charge, even less so to determine the structure of the 
representations. However, the finite-size transfer matrix is closely tied to the Temperley-Lieb algebra 
and allows us to push our examination much further. 

It was shown in [23] that the transfer matrix (actually its square) can be written in terms of spin 
operators (Pauli matrices) satisfying the Temperley-Lieb algebra TL n (/3) at fugacity /3 = 0. An imme¬ 
diate consequence is that the (spin) configuration space becomes a module over TL n (0). Partitioning 
the full space into sectors according to the value v of the total magnetisation gives rise to submodules, 
whose structures have been completely determined [23]. 

The presence of the algebra TL n (0) is key to assessing the conformal structure. Indeed, one can 
construct finite-size operators, as universal linear combinations of TL„(0) elements, which are believed 
to converge, in the scaling limit, to the Virasoro modes [47]. In this way, the TL n (0) modules, present 
at finite size, induce a conformal module structure in the infinite-dimensional space obtained in the 
scaling limit. We explicitly carry out this program in the present context and find, after a Jordan- 
Wigner transformation, operator modes Lk written in terms of fermion modes. We then verify that 
their commutation relations are those of the Virasoro algebra with central charge c = — 2. 

The scaling limit yields an infinite number of infinite-dimensional sectors E v , labelled by a quantum 
number v, each one being a Virasoro module. We determine the structures of these modules by 

2 The analysis made in [19] suggests that including trimers requires a conformal description with central charge c = 2. 
However, by extending the paradigm of dual interpretations, there might be more than one conformal description. 
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establishing an isomorphism between the E v and specific Feigin-Fuchs modules [48]. The result of 
this analysis is that each module E v , whose character is that of a single Verma module, contains an 
infinite number of (irreducible) composition factors, with a structure that depends on whether the 
scaling limit is taken over a sequence of odd or even lattice widths. In one case, the modules E v are 
direct sums of irreducible modules, while in the other, the composition factors form a single reducible 
yet indecomposable module (see Theorem 1 in Section 5.3). None of the modules exhibits non-trivial 
Jordan blocks for Lq. 

Finally, we examine the consistency of this conformal picture by linking it to the coefficients a p 
in the large n expansion of the energy levels. From (1.2), these coefficients are proportional to the 
eigenvalues of integrals of motion I r 2 P - i , which are, as discussed above, computable polynomials in 
the Virasoro modes. The action of the integrals of motion is completely determined by the conformal 
structure of the modules E v . Their eigenvalues hp-i can then be explicitly evaluated and used to 
compute the coefficients a p . In other words, the eigenstates of the transfer matrix, which converge 
to conformal states populating a certain level in a module E v , should be such that their associated 
coefficients a p , for each p, are related by (1.2) to the eigenvalues of Io^-i acting on that level. Comparing 
the so-obtained coefficients with their actual values computed from the lattice eigenvalues yields, in 
principle, an infinite number of verifications of the value of the central charge for each module E v . In 
practice, it is rather hard to perform them all, because the integrals of motion are not all explicitly 
known, or because the calculations are too cumbersome, as is the case for high levels. For the three 
lowest-lying levels in each module E v , and for the first two non-trivial integrals of motion, I 3 and I 5 , 
we find that there is complete agreement. 

These verifications constitute substantial evidence confirming that Lieb’s transfer matrix approach 
to the dimer model is consistent with a c = —2 conformal description. It is well known that different 
realisations of the Virasoro modes can be constructed in terms of free fermion modes, with different 
values of c. We take advantage of this freedom to construct, in terms of the same fermion modes as in 
the c = —2 description, an alternative realisation of the Virasoro algebra, this one with c = 1. Under 
this new conformal action, the full state space is partitioned into different subspaces E x , according 
to the values of a new quantum number x preserved by the Virasoro algebra. We determine the full 
structure of the subspaces E x as conformal modules and find a very different structure. Indeed, in all 
but one case, the E x are irreducible Verma modules. However, this alternative point of view turns out 
to be inconsistent with the transfer matrix description, in the sense that the eigenvalues of the integrals 
of motion do not reproduce the coefficients a p computed from the lattice. 

It is natural to speculate that a consistent conformal framework with c = 1 could nevertheless 
emerge from the use of a different transfer matrix, relying on different degrees of freedom. We hope to 
return to this question in the future. 

The paper is organised as follows. In Section 2, we recall the basics of the dimer model on a 
strip and its description in terms of Lieb’s transfer matrix. We give the details of its diagonalisation, 
which are useful in later sections. We also introduce the partitioning of the configuration space into 
sectors based on the total magnetisation and compute the corresponding partition functions. Section 3 
reviews aspects of the critical dense polymer model and a class of transfer tangles pertaining to the 
Temperley-Lieb algebra at /3 = 0. In Section 4, we first perform a detailed comparison of the spectra 
of the dimer and critical dense polymer models. We explain their similarities using a series of maps 
relating perfect matchings to configurations of the critical dense polymer model on a lattice tilted by 
45°. The integrability of the dimer model is also established in Section 4, where we exhibit a one- 
parameter family of matrices closely related to polymer transfer tangles and commuting with Lieb’s 
transfer matrix. Section 5 examines the way conformal invariance is realised by investigating in detail 
the two options c = — 2 and c = 1. These turn out to involve Feigin-Fuchs modules of the Virasoro 
algebra, whose definitions are reviewed in Appendix A. The proofs of the structures of the conformal 
modules are deferred until Appendices B and C. Section 5 also provides a detailed comparison of the 
integrals of motion and the lattice data, demonstrating the self-consistency of the c = — 2 description 
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and the inconsistency of c = 1 in the present context. Concluding remarks are presented in Section 6. 


2 Dimer model 

Here, we review definitions and results for the dimer model on the square lattice. The presentation 
follows [49] and [23]. 

2.1 Statistical model 

A dimer is a bridge that covers two neighbouring sites of a lattice. We define the model on the square 
lattice, where the sites are organised in an M x N rectangular grid, with M and N counting the rows 
and columns, respectively. Each site is labeled by a pair of integers (i,j) that gives its position in the 
grid, with 1 ^ i ^ N and 1 ^ j ^ M. A perfect matching is a configuration of Hp- dimers, where each 
site is covered exactly once. If MN is odd, the set of perfect matchings is empty. 

We fold the square lattice on a horizontal cylinder, meaning that the sites (*, 1) and ( i,M ), for 
1 ^ i ^ N, are neighbours and can be covered by the same dimer. In contrast, the sites (l,j) and 
( N,j ) are not neighbours and cannot be covered by the same dimer (unless N = 2). An example of a 
perfect matching on the 6x9 cylinder is given in Figure 1. 

The model depends on a free parameter a € R, which measures the preference for the horizontal 
dimers over the vertical ones. Perfect matchings are thus weighted by a h , where h is the number of 
horizontal dimers. The partition function is 

Z dim = J2 a h . (2.1) 

perfect 
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Figure 1: A perfect matching on the 6x9 cylinder and its corresponding spin configuration. 


2.2 Transfer matrix 

Lieb’s solution [8] of the dimer model uses a transfer matrix approach. The construction builds on 
a map from perfect matchings to two-dimensional spin configurations and ultimately expresses the 
cylinder partition function as the trace of the Af-th power of a spin-chain transfer matrix. 

A vertical edge is said to be occupied if a vertical dimer covers the sites at its ends. It is then 
assigned an up spin. Otherwise, it is unoccupied and assigned a down spin. This assignment of up 
and down spins to vertical edges maps a perfect matching to a two-dimensional spin configuration. An 
example is given in Figure 1. For a row of vertical edges (of length IV), we refer to its edge occupations 
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as a row configuration. Locally, the map sends a row configuration to an element of the canonical basis 
of (C 2 )® N , the vector space spanned by N ^-spins. For example, 




E (C 2 )® 4 . 


( 2 . 2 ) 


This map thus encodes the degrees of freedom of the vertical dimers. The occupation of a vertical edge 
is nevertheless fixed by the dimer attached to the site immediately below it. A site covered either by 
a horizontal dimer or by the top half of a vertical dimer always lies below an unoccupied edge. 

The canonical basis of (C 2 )®^ is made of the elements |s) = |siS2---sjv) with s* € {+,—} ~ 
{(o) > (?)}■ We denote by |d) the state with all spins pointing down: 


|d> = |-)■ 


(2.3) 


The Pauli matrices acting on the j-th component in (C 2 )®^ are 


with 


cr“ = I 2 <8> 


j - 1 


<g> I 2 _<8> <7 a ® I 2 <8> ■ ■ ■ <8> I 2 , 

N-j 


(a € {x,y, z, +,-}), 


(2.4) 


a 


X 







(2.5a) 


a = 





(2.5b) 


We now introduce Lieb’s transfer matrix T(a). For any |s), T(a) |s) is a linear combination of the 
states corresponding to the possible row configurations immediately above the row described by |s). 
T(a) is the product of two operators. The first one reverses all the spins, 


N 

Vi = U a J ( 2 - 6 ) 

3 = 1 

In terms of row configurations, acting with V\ on |s) produces a state (s') corresponding to a row 
configuration where vertical edges are occupied whenever possible, that is, unless the top half of a 
vertical dimer appears just below. For instance, 



(2.7) 


Of course, some perfect matchings also include horizontal dimers. If two adjacent spins in Is 7 ) are 
up, the state Is 7 '), where these two spins are flipped, corresponds to a row configuration of another 
perfect matching, where the two vertical dimers are replaced by a horizontal dimer that covers the two 
sites below. The second operation should thus replace |s 7 ) by a linear combination consistent with all 
possible insertions of horizontal dimers. This is implemented by acting on |s 7 ) with 


For instance, 


N -1 

^3 = H(I + a crj 
3 =1 


a 3+i) = 


N-l 

exp(]T 

3 = 1 


a a j a j+1 


( 2 . 8 ) 



+ Qi 





(2.9) 




























The result indeed consists of all the row configurations that can appear above ,Z 1 1 Z! . Although not 
illustrated in the previous example, V 3 can flip more than one pair of adjacent up spins if |s') permits 
it, thus allowing for more than one horizontal dimer in a given row of sites. It is worth noting that V 3 
also keeps track of the weights a. of the horizontal dimers. We arrive at the following expression for 
the transfer matrix: 

N—l N 

T(a) = V :i Vi = exp ( ^ a ajaj + ^j JJ <rj. ( 2 . 10 ) 

3 =1 3 =1 

For a € R, T(a) is real and symmetric, thus making it diagonalisable with real eigenvalues. T{a) 
is the appropriate transfer matrix for a strip domain and can be used to compute the partition function 
on the horizontal M X N cylinder (with vertical periodicity), 

z dim = Tr fM^ (2.11) 

The problem of computing Z dim is then reduced to the calculation of the eigenvalues of T(a), a matrix 
of dimension 2 N . 

We note that in general, 

[T(a),T(a')]^0, (2.12) 

implying that a family of operators that commute with T(a), and thus what we refer to as lattice 
integrability, is not naturally built-in. We will return to this question in Section 4.5. 


2.3 Diagonalisation of the transfer matrix 

In preparation for Sections 4.2 and 4.3, we restrict our attention to cylinders where the number of rows 
M is even. In this case, instead of T(a), we can study its square T 2 (a), which is conveniently written 
as 

N-l N -1 

T\a) = V :i Vj = exp ( ^ a ajaj +1 ^j exp ( a 0 + u+ +1 ). (2.13) 

3 = 1 i =1 

The calculation of the spectra and eigenstates of the transfer matrix was carried out in [ 8 ,49] using 
a Jordan-Wigner transformation. Here, we follow the same approach, but choose to first perform a 
unitary transformation that flips every odd spin, 


T(a) = UT(a)U~\ U= of. 

l^j^N 
j= 1 mod 2 


This yields 


LfJ LfJ 

T 2 (a) = exp ( a ( a Zj-i a 2j + (T 2 j a tj+ 1 )) ex P ( a ^ a 2 j-l a 2 j + a 2 j a 2 j+l)) . 
3 = 1 i =1 


where the convention cr ^, 1 = 0 is used. 
The hrst step in diagonalising T 2 (a) 


is to express it in terms of the operators 


Cj = (-1 y 


-1 


3- 1 

n 

k =1 


04 


a ■ 


c] = (- 1 ) 


3 - 1 


3 -1 

n < 

k =1 


U 


(1 < 3 < N), 


(2.14) 


(2.15) 


(2.16) 


which satisfy the usual (fermionic) anticommutation relations, 

{Cj , Cl } = S jtk , { Cj , c k \ = {Cj,C k } = 0 . (2.17) 
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With this transformation. 


LfJ LfJ 

T\ a ) = exp ( £ a (Ct^Cy + C\ j+1 C 2j )) exp ( £ a + C f 2J C 2J+I )). 


(2.18) 


J=1 


J = 1 


The next step is to perforin a Fourier expansion of C- and Cj, 


a = 


N 


E- 


sm | 1 r] C\ = 


V iV + 1 ^ “ ViV + 1/ 

9=1 


TV 


„ i / x sin - 7?'. 

3 \l N + 1 ^ VA1+ 1/ lq 

9=1 


E Sh 


This transformation is unitary, so the anticommutation rules are preserved: 

{V P , v\} = <W {V P , V g } = W P , v\} = 0. 


(2.19) 


( 2 . 20 ) 


In terms of these Fourier coefficients, T 2 (a) is expressed as a product of |_htJ mutually commuting 


blocks, 


L^J 


T 2 (a) = T q , T q = exp A+ exp A q , [T p ,T q \ = 0, 
9=1 


where 


,± ( nq 

At = a cos - 

q \N+ 1 


rtlvq - VN+l- q V N+1 - q ± (^jv+l-^9 “ ^TV+l-g 


( 2 . 21 ) 


( 2 . 22 ) 


The blocks T a can be diagonalised simultaneously. For N odd, Tjv+i = I is already diagonal. The 
H 2 
diagonalisation of the other blocks is obtained via the (unitary) transformation 

1 


Vq — /- (i^q "b 1) Cg + (Mg 1) C/v+i_g') ■ 

sj M +1) v ’ 

1 .) = ('7'... + 1) Cl + (P? - 1) Cf+l-g). 

\J2{Hq + 1) 


1?9 = 


(1 ^q^N), 


(2.23) 


where 


fi q = a cos 


\ . I n , 2 2 ( 

+ i /1 + cos- ( x 


JV+ 1 

Again, the commutation rules are preserved: 

(C P , 4} = <W {Cp, CJ = {Cj, cj} = 0 . 

Applying the transformation (2.23) to T q yields 

Tq = (nffiqtl-tN+l-qtN+l-q . 


(2.24) 

(2.25) 

(2.26) 


This T q is diagonal on the four-dimensional space spanned by |d), C||d), Cjv-i-i—cyI an< ^ Cjv+i- g Cg|d), 
with respective eigenvalues 1, ji 2 , 2 and 1. We note that for N odd, 

Cjv+l = V n+1 ) Cn+1 = 1 ?JV+1 • 

2 2 2 2 

The full transfer matrix squared is then obtained from the product of the blocks, 


(2.27) 




T 2 (a ) = (Mg )^ 9 <N+i-g^N+i - q . 

9=1 


(2.28) 
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The eigenvectors of T 2 (a ) take the form |v) = CgiC® • • • Cjvld)) and the corresponding eigenvalues, 
denoted by A, are obtained from (2.28) by replacing C<LC<?; by 1 for * = 1,... ,r and the other by 
0. After a change of variable, k = N + 1 — 2 q, they can be written in the following convenient form, 

A = (a sin 9 k + \J 1 + a 2 sin 2 Q^j \ S k = u k + r k , 0 k = 0 ^' + ^ (2-29) 

k=N—l mod 2 


where the freedom in the presence of the excitations Q is now encoded in the quantum numbers 
v k , r k € {0,1} defined by 

v k\ v ) = C iv+i+fc C iv+i+fc I v )) 

2 2 (1 < IV-1, A; = JV - 1 mod 2). (2.30) 

r fc| v ) = (1 — CiV+l-fcCiV+l-fc)l V )l 


For N even, there are 2 N choices of vectors {y,r). For N odd, each of the 2 N 1 vectors {y,r) appears 

twice due to the block Tn+i = I, and the total number of eigenvalues is again 2 N , as required. 

2 

For a > 0, the largest eigenvalue 3 of T 2 (a), denoted by Ao, corresponds to u k = Tt = 0 for all k. 
It is non-degenerate for N even, but appears twice for N odd. From (2.29), the energies of the ground 
state and of the excited states are found to be 


ZT’dim 

^0 




ZT'dim 

^0 


N—l 

— |logAo= — ^2 arcsinh(a sin 8 k ), 
k= l 

k=N—l mod 2 


-iiog 


A_ 

Ao 


JV—1 


(v k + Tfc) arcsinh(a sin 0 k ). 


k =1 

k=N— 1 mod 2 


(2.31a) 


(2.31b) 


2.4 Total magnetisation and variation index 


The transfer matrix and its square, respectively, anticommutes and commutes with the total magneti¬ 
sation S z , 

N 

{T(a),S z } = 0, [T 2 (a),S z } = 0, £* = £!>?. (2.32) 

3=1 

This was first observed at the level of T(a) and T 2 (a ) in [49], where the operator 

V = U~ l S z U 


was called the variation index. Each of the eigenspaces 


E° n = (C 2 )®^' 




with dimension 


dim A’y = 


r r N N — l N—l n -\ 

u t 2 ’ 2 ’■■■’> 2 I 2 1’ 

N 


N 


— V 


(2.33) 

(2.34) 

(2.35) 


subspaces for v = ±4y- is seen, at the level of the dimers, as coming from the (separate) invariance 


is thus invariant under the action of T 2 (a). In particular, the existence of the two one-dimensional 
subspaces for v = ±y is seen, at the leve 
under T 2 (a) of the two row configurations 

J o o o o 

o — • - 0 —!• - 


and 


o o o o o 

O - • — o - .| o 


(2.36) 


3 Because A (—a) = A(a) 1 , the maximal eigenvalue of T 2 [a) for a < 0 instead corresponds to Uk = Tk = 1 for all k. 
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(2.37) 


The total magnetisation can equivalently be expressed in terms of the fermions as 

N N N 

^ = E ( c ) c j - P) = E fak - P) = E (cjc, - P)- 

j =1 9=1 9=1 

The one-dimensional spaces mentioned above are spanned by |-- ■ ■ —) for v = — ^ and | + + • • • +) 

for v = y. 

Only a subset of the eigenvalues of T 2 (a) appear in each sector E V N . An eigenvalue is thus 
v -admissible if its eigenvector has eigenvalue v under S z : C 91 C 92 • • • C|r|d) is in the sector E V N with 
v = r — y- Equivalently, in terms of the quantum numbers and Tk, an eigenvalue A as given in 


(2.29) is u-admissible if 


n-i r 

E k - r fc) = { ,1 
\ V+ 2 


or v — 


N even, 
N odd. 


(2.38) 


k=N—l mod 2 


In particular, the two degenerate eigenvectors of T q discussed below (2.26), |d) and CjCjv+iE^’ belong 
to different magnetisation sectors. 

2.5 Partition functions 

In the scaling limit, M and N are taken very large with the ratio M/N converging to a finite non-zero 
number 5. The finite-size corrections (1.1) are extracted from (2.29) using the Euler-Maclaurin formula. 
The bulk and boundary free energies are given by 4 




1 


arcsinh(a sin t) d t, /bd y = - arcsinh a, 


(2.39) 


/bulk — " 

whereas the speed of sound and the effective central charge are 

f -2 IV odd, 

0 = «> c eff = < , AT 

I 1 J\ even. 

The full partition function can be split into contributions from the different E V N sectors, 

z dim = ^ aM/2 > Z dim = ^zp im . 

u-admissible A v 

The conformal partition functions Z dim (q) and Z dim (q ) are defined as 

e M(Ar+l)/ bulk +M/ bdy ^dim M,N^oo^ 2 dim (q) e M(lV-H)/ bu ik+M/ bdy ^dim M,N^o ^ 42) 

It was shown in [49] that the Z dira (q) are equal to Virasoro Verma module characters, 


(2.40) 


(2.41) 


Z d ™{q) = 


q 


,u 2 /2 


r){q) 


q = e 


—anS 


(2.43) 


for both parities of N. The full conformal partition function Z dim (q) is obtained by summing Z d ' m (q) 
over v € Z for N even and v £ Z + ^ for N odd. This yields 

03 (q) 


Z dha (q) = '£Z? m (q) 


= 


v(q) 

02(g) 

. v(q) 


N even, 
IV odd, 


(2.44) 


4 We note that for the dimer model, Hnite-size corrections of the form (1.1), with only odd powers of / in the sum, 
require that we set n = IV + 1. This explains why the boundary free energies given here and in [49] are different. 
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where 9j(q) is the j-th standard Jacobi theta function and r](q) is the Dedekind eta function: 


1 

0 3 (q) = ^2qj 2 / 2 , v (q) = q2i J|(l - q J ). (2.45) 

jez jez j =1 

3 Critical dense polymers 

This section reviews definitions and results for critical dense polymers [43]. 


3.1 Statistical model 


The model of critical dense polymers is defined in terms of nonlocal observables called loop segments. 
We consider a lattice constructed as an m x n array of square faces, with m even. Periodicity is 
imposed in the vertical direction, so the lattice is folded on a horizontal cylinder. Each of the mn faces 
is decorated by either Q or The midpoint of a face edge is called a node and this is where the 
loop segments are attached and connected to those of neighbouring faces, thereby forming longer loop 
segments in the form of curves on the surface of the cylinder. At both ends of the cylinder, we draw 
additional loop segments in the form of semi-circles that connect each odd row with the even one just 
above it. The resulting diagram is a loop configuration a. Figure 2 presents an explicit example. 



Figure 2: A configuration of critical dense polymers on the 6x9 cylinder. 


Every loop segment is part of a (closed) loop which is either contractible or non-contractible. 
For the model of critical dense polymers, the fugacity ft of contractible loops is zero, meaning that a 
configuration a with one or more such loops is assigned the weight W(a) = 0. Other configurations, 
those without contractible loops and with n 7 non-contractible loops, are weighted by 

W{a) = 7 n ^pf ' m pf :S) . (3.1) 

Here, 7 € M is the fugacity of the non-contractible loops, and p\,p 2 € K are local weights associated 
with the two types of 1 x 1 diagrams. More generally, p\ and P 2 can be position-dependent and allowed 
to take different values on the various faces. Finally, the partition function is the sum of the weights 
over the set of configurations: 

Z cdp = J2 W ( (J )- (3.2) 

(7 


3.2 Temperley-Lieb algebra 

The Temperley-Lieb algebra [50,51] on n sites, TL n (/3), appears naturally in the model of critical dense 
polymers because its rules for products of connectivities mathematically encode the nonlocal rules that 
define the loop model. The parameter (3 € C is the fugacity of contractible loops. For critical dense 
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polymers, contractible loops are disallowed, so (5 = 0. The representation theory of TL n (/3) was studied 
by Martin [52], Goodman and Wenzl [53], and Westbury [54], see also the recent review article by 
Ridout and Saint-Aubin [55]. 

The elementary objects spanning the algebra TL n (/3) are connectivities. Let us draw a rectangle 
with n distinct nodes on the top edge and n more on the bottom one. A connectivity is a pairwise 
connection of these nodes by non-intersecting loop segments, all drawn inside the rectangle. For 
instance, 


is a connectivity in TI_7(/3). Two connectivities are considered equal if their nodes have the same 
connections. The number of distinct connectivities is given by 


(3.3) 



dimTL n {fi) = —- f 2n ^j . 

n + 1 V n J 


(3.4) 


Addition is commutative and the linear combinations of connectivities over C are called tangles. The 
algebra TL n (/3) is the vector space spanned by the connectivities, endowed with the following product 
rule. Let a\ and <22 be two connectivities in TL n (/3). The product a\a ,2 is obtained by drawing a\ under 
02, superposing the top edge of a\ and the bottom edge of in such a way that the n nodes of the two 
edges match up. The edge common to aq and 02 is then removed. This produces a bigger rectangle 
where the nodes from the top and bottom edges are connected, thus giving rise to a new connectivity 
03. If closed loops are formed, they are erased and each is replaced by a multiplicative factor of (3. For 
example, 


a 1O2 — 



= /3 2 







= /3 V 


(3.5) 


The algebra TL n (/3) is alternatively defined in terms of generators, 


TL n (P) = (I, ej] j = 1 ,... ,n — l), 

satisfying the relations 

19 = 91 = 9, ej = (3e r 

where g is any of the generators. 



123 n 




1 j j +1 n 


(3.6) 


&i&j — Cj Cj (|% jj 1), (3.7) 


3.3 Standard representations 

Computing physical quantities in the loop model typically requires working with representations of the 
Temperley-Lieb algebra. A family of representations that appears naturally, the standard representa¬ 
tions, is based on the construction of link states. Let there be n distinct nodes on a horizontal line. 
A link state is a diagram where each node is occupied by a loop segment living above the horizontal 
line. These non-intersecting loop segments either connect nodes pairwise or are vertical loop segments, 
called defects, that cannot be overarched. For obvious reasons, the number d of defects is constrained 
to have the same parity as n. The linear span of link states on n nodes with d defects is denoted by 
V((. For example, 


The action of a connectivity a on a link state w G V(( closely resembles the rule given for the 
product of two connectivities. To compute aw, one draws w above a, erases the common horizontal 
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line segment, reads the new link state from the bottom n nodes and adds a factor of /3 for each 
contractible loop created (and subsequently erased) in the process. The standard action has an extra 
rule: If the number of defects has decreased, the result is set to zero. For instance, 



For each 0 ^ d ^ n with d = 
dimension is 



II xTToh 

fd 1 r\ 1 , 

KJ) K\JJ 

frr^\ r\ 


n mod 2, this defines a representation of TL n (/3), denoted by pd, 


(3.9) 

whose 


dim pd = dim \/ d 


n 

n—d 

2 


n 

n—d—2 


(3.10) 


3.4 Gluing operator 

From here onwards, we focus exclusively on the case fd = 0. Forgetting for the moment about the vertical 
periodicity imposed on the lattice, configurations of the loop model can be viewed as connectivities of 
TL n (fd = 0). The example in Figure 2 corresponds to 


a = 



(3.11) 


By splitting the loop configuration into the three 2x9 configurations corresponding to the connectivities 

(3.12) 


a i = 



j a 2 ~ 



> «3 — 

f r\ r\ 


the connectivity a can be expressed as the product a = 010203. 

The partition function Z cdp can be computed using the linear operator 


T : TL n (0) -7 C. (3.13) 

For each connectivity a, dF{a) outputs the product of the fugacities of the loops produced by folding 
a around the cylinder and gluing together its top and bottom edges. We therefore call T the gluing 
operator. If the folding produces contractible loops, then F(a) = 0. If only non-contractible loops are 
formed, then F(a) = with 7 and n 7 defined below (3.1). For example, for (3.11), F(a) = 7 s . 

One known realisation of T is in terms of traces in the standard representations [56,57]: 

n 

Ha)= E U d Q)Trp d (a), (3.14) 

d=0 

d=n mod 2 

where U^ix) is the k -th Chebyshev polynomial of the second kind. The linearity of the realisation 
(3.14) of T is inherited from the same property of the matrix trace. 


3.5 Transfer tangle 

It is possible to encode all the connectivities needed for the calculation of Z cdp in a single TL n (0) 
tangle, the double-row transfer tangle It is defined as 


D(u,0 = 


1 r 

u— £i 

U — £2 
r 



U — ^n 

r 

sin2u 

1 

^_ 

^_ 



U + £n 

^_ 


(3.15) 
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where each square tile is called a face operator and is formally given by 


u 

= cosit 

1 _ 



A 


+ sin u 


(3.16) 


The small quarter-circle in the lower-left corner of is a marker indicating the orientation of the 
diagrams in the decomposition. A face operator with its marker in another corner is obtained by the 
corresponding rotation of the diagrams on the right-hand side of (3.16). The parameter it € C is the 
spectral parameter while £i,£ 2, • • • ,£n € C are inhomogeneities. We often use the shorthand notation 
£ = (Ci) £ 2 , • • • j £n)j as in the left-hand side of (3.15). We note that the normalisation of D(u, £) ensures 
that the Oth order term in its Taylor series expansion in u is non-zero. For £ = 0, the resulting term is 
a multiple of the identity connectivity. 

As a TL n (0)-tangle, D(u, £) is defined in a diagrammatic way. Each face operator has two contri¬ 
butions from which one obtains a linear combination of 4” diagrams. Each of these takes the form of a 
2 x n loop configuration from which a connectivity can be read off. Some of these configurations contain 
contractible loops and are discarded because j3 = 0. D(u,£) is defined to be the linear combination 
of the remaining ones, each weighted by the corresponding powers of p\ and p 2 - here parametrised 
by the position-dependent functions cos(it ± £j) and sin(tt ± £j). Two transfer tangles with identical 
homogeneities commute [58,59], 


[D(u,0, D (v,£)] = 0, u,v 1 E C. (3.17) 

The parameter u in D(u, £) thus explores a family of commuting matrices — accordingly, it is usually 
referred to as spectral. In stark contrast, -D(u,£) and D(u , £') do not commute in general if £ / £', and 
the inhonrogeneities are therefore not spectral parameters. The multi-parameter definition of Z}(u,£) 
provides a considerable freedom in choosing the position-dependent weights p\ and P 2 in (3.1) while 
ensuring that the model is integrable. 

The tangle (.D(ii, £)) m ^ 2 is a linear combination of connectivities, each weighted by appropriate 
powers of local weights such as p\ and po in (3.1). They include those produced from the loop con¬ 
figurations contributing to Z cdp , but also some which are mapped to zero by the gluing operator. 
Applying T to this linear combination removes the unwanted connectivities and correctly assigns y ” 7 
to the surviving ones, so that 

Z cdp = T(D(u^) m/2 )- (3.18) 

From (3.14), the calculation of Z cdp amounts to finding the eigenvalues of the transfer matrices 
Pd(D(u, £)) for each d. 


3.6 Spectra 

The technique used to find closed forms for the eigenvalues of pd(D(u, £)) mimics a method originally 
used to obtain the spectrum of the Ising model transfer matrix [60,61]. The idea is to show that two 
transfer matrices, with their spectral parameters related in a specific way, are inverses to one another, 
up to a scalar function of these spectral parameters. The eigenvalues are then computed from this 
scalar function. 

For critical dense polymers, an inversion identity indeed exists. In fact, it is satisfied at the level 
of the algebra TL n (0), in the sense that the product of two specifically related transfer tangles gives a 
scalar multiple of the identity connectivity. This was first shown in [43] for £)(«,£) for £ = 0 and later 
generalised to other geometries and boundary conditions [62-64]. The ideas used for D(u, £) at £ = 0 
easily extend to the general case where £ is free. The result is 


£>(u, £)£>(« + §,£) = / 


f fl(u,Q - f2(u,Q 

\ cos 2 u — sin 2 u 


2 


(3.19a) 
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with 


fi(u,0 = JJcos(u-£j)cos(u+ £,■), f 2 {u,€) = JJsin(u-^)sin(u + ^). (3.19b) 

3=1 3 =1 

In any representation, each (parameter-dependent) eigenvalue Eig(Z?(u, £)) of D(u , £) then satisfies 
the relation 

Eig(Z?(u,£))Eig(r>(u+ §,£)) = . (3.20) 

V cos u — sm u J 

Factorising the right-hand side of (3.20) is not possible in the general case, but it can be done for 
special choices of £. For £ = 0, this factorisation yields the following expressions for the eigenvalues, 


L^J 


Eig^M)) = n 

3 = 1 


1 + €j sin 2 u sin tj'j ^1 + pj sin 2 u sin . 




pr_ 

n 


2 n 


n even, 
n odd, 


(3.21) 


whose forms are simplified here compared to those in [43] . Different eigenvalues correspond to different 
choices of ej, p 3 G {+1,-1}. With v\. = and r k = l ~J ? k , these can be recast into 


n—2 

Eig(D(it,0)) = n ( 1 — sin 2 2 u sin 2 q2j 

k =l 

k=n mod 2 


1 + sin 2 u sin q k 
1 — sin 2 u sin q k 


5k 


v k + T~k, 


irk 

Qk = 2 ^’ (3 ' 22) 


where the eigenvalues are now specified by a selection of v k ,T k £ {0,1}. Eigenvalues of D(u, 0) are 
denoted by A. For sin2n > 0, the largest one, Ao, is obtained by setting v k = T k = 0 for all k, for both 
n even and odd.° Because D(u, 0) covers two rows of n tiles each, we include, by convention, a factor 
of | in the definition of the energies: 


n—2 


TP cdp 
^0 


E cdp _ E cdp 


- k log Ao = - Y log(l + Sin2ri sin^fc), 


k=1 

k=n mod 2 
n— 2 


h°e(y)= S &log VT 


1 + sin 2 u sin q k 


k =i 
k=n mod 2 


— sin 2 u sin q k 


(3.23a) 


(3.23b) 


In a given standard representation pd , there are restrictions on the possible values that the quantum 
numbers v k and r k can take (otherwise, the number of eigenvalues would exceed the dimension of 
Pd{D(u, 0))). One therefore needs to state which pairs of vectors (z/, r) are retained to form the 
eigenvalues in a given standard representation. The appropriate admissibility criterion was conjectured 
in [43] and later proven in [65] , and is more complicated than the one for the dimer model. This criterion 
states that a pair (y. r) is d-admissible, that is, defines through (3.22) an eigenvalue of pd{D(u, 0)), if 
its components satisfy the two following constraints, 


Y, (vk-Tk)> 0 (AT N), 

k^K 

k=n mod 2 


n—2 


Y ("k~T k ) 

k= 1 

k=n mod 2 



d-2 
2 


d-1 
2 


n even, 
n odd. 


(3.24) 


Although the values v k , r k are severely constrained, all the possible eigenvalues (3.22) appear in at 
least one standard representation, that is, every 5 with 5 k G {0,1, 2} corresponds to a d-admissible pair 
(v, t) for at least one value of d. 

J For sin2u < 0, the maximal eigenvalue instead corresponds to Vk = Tk = 1 for all k. 
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3.7 Partition functions 


Equipped with the closed-form expressions for the eigenvalues and their degeneracy criteria, one can 
calculate the cylinder partition functions restricted to each p d as well as the full partition function: 

n 

Z c d dp = Tr(p d (D(u,0) m / 2 )) = mult a \ m / 2 , Z cdp = ]T U d Q) Z c / P . (3.25) 

d-admissible A d =0 

d=n mod 2 


We now study the behaviour of the partition function as m, n —> oo and ^ —>• 5, with 5 € R. The 
continuum scaling limit will be discussed further in Section 5.1. By performing an Euler-Maclaurin 
expansion of E^ dp in (3.23a), the bulk and boundary free energies are found to be 


1 f 71 "/ 2 1 

/bulk =-/ log (1 + sin 2 u sin t)dt, f h dy = - log (l + sin 2 u ), 

71 Jo 1 

while the speed of sound and effective central charge are 


i? = sin 2 u, 


c e ff = 


—2 neve n, 
1 nodd. 


(3.26) 


(3.27) 


The conformal partition functions are obtained by renormalising the partition functions to remove 
the non-universal terms and by taking the scaling limit: 


e mn/ bu lk+m/bdy ^dp m ’"^°°) Z^ dp (g), e mrl /bulk+"l/bdy Z cd P Z cdp (q). 


(3.28) 


For each d, Z d dp (q) produces the so-called Kac character Xi,d+i(q) [43,59], 

(d- 1) 2 

Zf V {q) = Xi, d+ i(q) = ^-(1 - q d+1 ), q = e~^ 2u . (3.29) 

From (3.25), this gives an explicit expression for the full conformal partition function Z cdp (q), for all 
7 . Because U d { 1) = d + 1, the result simplifies significantly for 7 = 2 : 


d 


' 2 d 2 (q) 

v(q) 
203(g) 
. v(q) 


n even, 


n odd. 


(3.30) 


The similarity between (2.44) and (3.30) is striking and is discussed in Section 4. 


4 Dimers, loops and integrability 

One of our goals is to study the conformal properties of the dimer model. In this section, we initiate 
this investigation by establishing a relationship between the dimer model and the critical dense polymer 
model, as finite-size lattice models. 


4.1 Comparison of spectra 

A priori, the models of dimers and critical dense polymers are quite different. The objects defining 
the former, the dimers, cover only two neighbouring sites and are in this sense local objects, while the 
polymer model is based on nonlocal entities, the loop segments, that can extend over long distances 
but not form contractible loops. Each model can be solved by analytically computing the eigenvalues 
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of a transfer matrix, although the techniques employed to find their closed forms are different. Another 
distinguishing feature is that the loop transfer tangles D(u, 0) and their matrix representatives form a 
commuting family, while this is not the case for the dimer transfer matrices T(a). 

Nevertheless, one sees a formal but striking similarity in the finite-size spectra (2.29) and (3.22) 
of T 2 (a) and D(u,0) provided the system sizes are related as 

n = N + 1. (4.1) 


In what follows, we will assume this relation. As observed in [49] , the similarity of the spectra strength¬ 
ens when we look at the finite-size corrections of the energies. In both models, the ground-state energy 
is of the form 


n —2 

E o = - w(g fe ), 

k=1 



(4.2) 


k=n mod 2 


but with model-dependent functions u(t), 


w dim (f) = arcsinh(asinf), a; cdp (t) = log (1 + (sint), 


(4.3) 


each depending on an extra parameter, a or £, where 

£ = sin 2 u. 


(4.4) 


In both models, the A expansion of Eq is computed using the Euler-MacLaurin formula and 
expressed in terms of the series expansion of the corresponding 


At, 




P= 1 


p\ 


(4.5) 


The ground-state energy Eq then has an expansion of the form (1.1) in the variable n, with 

1 


/bulk =-/ u(t)dt, 

7T Jo 


(o) _ ’r 2p lB 2 P (r) 

a p ~~ 


(2 P)l 


fbdy ~ 2 ^ 2 )’ 

A 2p -i, r = § mod 1, r € {0, |}, 


where B 2 p(z) are Bernoulli polynomials. Useful values are 


1 


1 


1 


B 2(0) = B 4 (0) = ~ B 6 (0) = ^, ■■■ and B 2p (i) = (^^—f - l)B 2p (0), 


(4.6a) 

(4.6b) 

(4.7) 


where the B 2p (0) are the Bernoulli numbers. It is noted that, although the series expansion of ui(t) 
may involve both even and odd powers, the large n expansion of Eq depends only on the coefficients 
of the odd powers. For p = 1, setting d = \\ and comparing with (1.1) yields the value of the effective 
central charge, 

—2 neven, 

1 nodd. 

The analysis of the energy gaps is similar. In both models, they are given by 

n —2 

A E = E - E 0 = ^ \ {v k + r k ) [u(q k ) - uj(-q k )]. 


c e ff = —12B 2 (r) = 


(4.8) 


(4.9) 


k=1 

k=n mod 2 


The odd combination u(t) — uj(—t ) in this expression implies that the large n expansion of the energy 
gaps, as in the case of Eq, only depends on the odd-indexed coefficients A 2p -1 ■ In fact, for the dimer 
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model, ui(t) is an odd function, so the antisymmetric combination " (/| —— coincides with u(t). As 
the finite energy excitations involve only a finite number of non-zero n k and rfc, it is not difficult to see 
that the A expansion for the energy gaps yields 


cip 



/7T\ 2 P -1 Agp-l 
V 2 / (2p — 1)! 


^ 1 {U k + T k ). 

k^l 

k=n mod 2 


(4.10) 


Putting it all together, we find that, in both models, the energy levels are unequivocally specified 
in terms of a set of admissible numbers {nki^k} and an appropriate function, cj dim or w cdp . Albeit 
different, the finite-size corrections in the two models are structurally identical to all orders in terms 
of n and the corresponding function uj. Explicitly, the expansion coefficients read 


7r 


2p-l 


dn — 


(2p-l)! 


A2p-1 


B 2 p(r 

2 p 


+ 


2 2 P~ x 


£ k^in + Tk) 


k=n mod 2 


(4.11) 


and the only model dependence is in the factors A 2 p-i • These factors are odd polynomials of degree 
2 p — 1 in their respective parameter, but are distinct, as the first few values demonstrate: 


Xf m (a)=a, A^ im (a) = ~a(a 2 + 1), = a(9a 4 + 10a 2 + 1), (4.12a) 

Af p (C) = C, A c 3 dp (C) = C(2C 2 - 1), Af p (C) = C(24C 4 - 20C 2 + 1). (4.12b) 


We see that the finite-size correction coefficients a p are not universal. However, the ratios a p /\ 2 P ~i, or 
equivalently the ratios a p /a^\ are believed to be universal. The universal nature of the similar ratios 
for the Ising model and critical dense polymers was discussed in [27] and [66], respectively. 

Up until now, the analysis of the finite-size corrections has been performed without any assump¬ 
tions regarding the central charge of the underlying conformal field theory. In the remainder of this 
section, we will work under the hypothesis that the dimer model, like critical dense polymers, can be 
described by a conformal held theory with c = — 2. We will return to this issue in Section 5. 

With this assumption, it follows that the set of distinct conformal weights in the scaling limit, 
computed from the coefficients ai, 

A = C 2 ^° ff + \ k ( u k + r k ), (4.13) 

fc^l 

k=n mod 2 


is exactly the same in the two models. The two conformal spectra, however, are not identical, because 
the eigenvalue degeneracies are different. For instance, for fixed d and v, two pairs (nk,Tk) and (n' k ,T' k ) 
may be such that v k + T k = u' k + t! for all k, with one being both d- and u-admissible while the other 
is only u-admissible. In fact, the admissibility criteria (2.38) and (3.24) show that the vectors (u, r) 
admissible in V d are a subset of those admissible in E^_ 1 , v = 

This can also be seen from their conformal partition functions, when the two models are defined 
on a cylinder of length n = N + 1 and perimeter m = M. The conformal partition functions of both 
models can be written in terms of the c = —2 irreducible Virasoro characters 


/ -i T 

ch r,s(9) = 9 J 3r +Ar, ‘ 9 


v(q ) 

where r, s are Kac labels and the conformal weight is 

(2 r — s) 2 — 1 


(r G Z^o, s = 1,2), 


A r,s — 


8 


(4.14) 


(4.15) 
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For the two models considered here, the explicit forms of the cylinder partition functions strongly 
depend on the parity of n. They read, for n even, 

oo 

Zv im {q) = Y ch FI+r+i,i(^) ( v € Z + i)> zf p (q) = ch d A (q) + chd+i^q) ( d € 2Z), (4.16a) 

r =0 

whereas, for n odd, 

OO 

= £ ch M+ 2 r+ 1 , 2 (g) («€Z), ^ dp (g) = clm+i ,(g) (d€2Z + l), (4.16b) 

r=0 2 

where we note that cho.i (q) = 0 . 

Sections 4.2 to 4.4 provide a deeper understanding of the similarities between the spectra, by 
building a direct connection between the two models. We also emphasise at this point that the conformal 
partition functions are mere generating functions for the conformal spectra; they contain very little 
information about the structure of the Virasoro representations making up the spectra. We will come 
back to this in Section 5.3. 

4.2 Dimers and webs 

Starting with Temperley’s correspondence [24] between dimers and spanning trees, we now establish a 
direct connection between the dimer model and critical dense polymers. While the authors of [42] have 
used the transfer matrix formulation of the dimer model to probe combinatorial aspects of spanning 
trees and spanning webs (see below), here we exploit Temperley’s correspondence to understand the 
similarity of spectra described in Section 4.1 and ultimately to study the conformal properties of the 
dimer model. Along the way, it will yield a solution of the dimer model entirely in terms of a loop 
model. The organisational diagram of maps used in this section and the next is displayed in Figure 3. 


perfect 

matchings 


oriented 

webs 


unoriented 

webs 


dense polymer 
configurations 


Temperley-Lieb 

connectivities 


Figure 3: The organisational diagram of maps from perfect matchings to TL„(0) connectivities. 

The first arrow in this diagram is Temperley’s correspondence first discussed in [24], The construc¬ 
tion is applicable to lattices with a variety of boundary conditions, but here we focus on the horizontal 
M x N cylinder defined in Section 2.1. 

Starting with the MxN square lattice, one can define a sublattice made of the sites (i, j) with i = 
j mod 2. In any perfect matching, each dimer covers one site of the sublattice. Instead of characterising 
dimers with regard to their horizontal or vertical orientations (as in Figure 1), we now distinguish them 
according to the parity of i (and j ) of the site (i,j) of the sublattice that it covers. In the example 
of Figure 4 (a), dimers on odd and even sublattice sites are, respectively, coloured in green and red. 
Requiring that the sublattice be periodic in the vertical direction constrains M to be even, thus 
justifying this choice already made in Section 2. 

For a given perfect matching, Temperley’s correspondence is constructed as follows: For each site 
(i,j) of the sublattice, one draws an arrow of the appropriate colour, starting from (i,j) and pointing 
to — (i+2,j), (i,j — 2) or (i,j + 2), following the orientation of the dimer occupying Here, 

j is taken modulo M, so that arrows starting at (i, 1) and (i, 2) can, respectively, point to ( i , M — 1) 
and ( i,M ), and vice versa. For a site (i-j) of the sublattice with i = 2 or JV — 1, it may happen 
that the arrow produced by this prescription points to a site outside of the lattice, in position ( 0 ,j) or 
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(N + 1 ,j). These external sites are referred to as roots and drawn using white circles in Figure 4 (b). 
For N odd, the roots on both sides of the cylinder are reached by red arrows, whereas for N even, the 
roots on the left and right are reached by red and green arrows, respectively. 

In the diagram that emerges, called an oriented web , the sites of the sublattice are connected 
by arrows that either (i) flow to a root, producing a rooted tree , (ii) are part of a unicycle, that is, 
a chain of arrows that wraps around the cylinder, forming a non-contractible cycle, or (iii) flow to 
a unicycle. Red and green structures are intertwined and space-filling [42], but non-intersecting. For 
obvious reasons, the number L of unicycles of an oriented web is constrained to have the same parity as 
N. Other important features are that there are no free ends with arrows pointing towards them, except 
at the roots, and that there are no contractible cycles of arrows. Indeed, such formations would cover 
an odd number of sites of the original M x N dimer lattice and thus never appear from Temperley’s 
correspondence, simply because any chain of dimers covers an even number of sites. 

Starting from an oriented web, it is not hard to reconstruct the unique perfect matching in the pre¬ 
image. Temperley’s correspondence is therefore bijective. Moreover, the number of horizontal dimers 
h is equal to the number of horizontal arrows in the oriented web. Temperley’s correspondence thus 
provides an alternative way to calculate Z dim : 

Z dim = J2 a h . (4.17) 

oriented 

webs 

By removing the orientations of the arrows, one obtains an unoriented web, composed of line 
segments called branches, all of equal length (twice the lattice spacing). This is the second map in 
the organisational diagram of Figure 3. The concepts of trees and unicycles carry over to unoriented 
webs in the obvious way. In the example of Figure 4, the unoriented web appears on panel (c), where 
a loop configuration relevant to the discussion of Section 4.3 has also been superimposed. The map 
from oriented webs to unoriented webs is clearly not one-to-one, but it is surjective. Indeed, given an 
unoriented web with L ^ 1 unicycles, there are more than one oriented web in the pre-image. Other 
oriented webs mapped to the same unoriented one are obtained by reversing the orientation of all 
arrows in one or more of the unicycles. At the level of the dimers, this is equivalent to shifting by one 
position the dimers in the selected unicycles, while leaving the other dimers unchanged. We conclude 
that an unoriented web with no unicycles has a unique pre-image, and, because each unicycle can be 
oriented in two ways, an unoriented web with L unicycles has a pre-image of dimension 2 L . 

The partition function of the dimer model can then be expressed in terms of sums of unoriented 
webs. The number h of horizontal dimers is now given by the number of horizontal branches. Following 
the previous discussion, to reproduce the correct counting of the oriented webs in (4.17), each unoriented 
web must be weighted by an extra factor of 2 L . This yields 

Z dim = J2 a h 2 l . (4.18) 

unoriented 

webs 


4.3 Webs and loops 

In the previous section, perfect matchings were mapped to unoriented webs made of branches that 
never form contractible cycles. This is reminiscent of the rule imposed on loop configurations in the 
critical dense polymer model, where contractible loops are disallowed. We now establish a direct link 
between unoriented webs and loop configurations. This is the third map in Figure 3. 

The first step is to superimpose, on the unoriented web, a new square grid, tilted at a 45° angle, 
that draws edges between each site (i, j) of the even sublattice and the site’s four neighbours (*±1, j± 1). 
Extra edges are added to connect each site (l,j) with j odd to the two roots in positions (0, j ± 1), 
and similarly for ( N,j ) and (N + 1, j ± 1) with j = N mod 2, while taking the vertical periodicity into 
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Figure 4: The results of mapping a given perfect matching to an oriented web, an unoriented web and 
a loop configuration. 


account. The result is a grid of Hp- diamond-shaped tiles covering the cylinder. Each of these tiles is 
crossed by a branch that connects two opposite corners. To establish the relation with the loop model, 
we decorate each tile with two quarter circles using the rule 



(4.19) 


for branches of either colour. By adding the boundary triangles 



(4.20) 


on the left and right edges of the tilted grid, we obtain a diagram where the quarter circles form 
loop segments that trace the contours of the structures of the unoriented web. This is illustrated in 
Figure 4 (c). Because contractible cycles are disallowed and each tree is attached to a root, the loop 
segments never form contractible loops. Instead, they produce a loop configuration of critical dense 
polymers on the cylinder, but on a tilted lattice compared to the one described in Section 3. 

Starting from the left, the unicycles in any unoriented web appear alternately in green and red, 
separated by single non-contractible loops. For any configuration with at least one unicycle, the loop 
configuration it maps to contains two non-contractible loops that act as outer contours for the leftmost 
and rightmost unicycles. The number of non-contractible loops is thus n 7 = L + 1. This also applies 
in the case L = 0 (and N even) where there are no unicycles and only a single non-contractible loop 
separating the left (red) and right (green) trees. 

This map from unoriented webs to loop configurations is readily seen to be bijective, and as the 
branches do not provide any information not already encoded in the contours, they are now erased. 
The partition function for dimers can thus be computed using weights W (u) for loop configurations a 
on the tilted lattice: W(a) = 0 if a contains contractible loops, and 

W(a) = 7 (7 = 2), (4.21) 

otherwise. Because #(<$>) + #(<$>) = W(<j) only depends on the ratio P 1 /P 2 up to an overall 

factor. For fixed a, the values of p\ and p 2 for the dimer model are given by 

a 1 , 

Pi = -, P'2 = (4.22) 

K K, 


for some value of the otherwise free parameter n. From (4.18), we find 

Z dim = k MN/2 ^ 2 ^^pf^pf^ = 4 K MN / 2 Z cdp 


7=2 


(4.23) 
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where Z cdp is the critical dense polymer partition function on the tilted lattice. 6 7 

From Sections 3.4 and 3.5, Z dim can be computed using a transfer tangle in TL n (0) and its matrix 
representatives in the standard representations pd . Omitting the folding on the cylinder, the loop 
configuration produced from each unoriented web can be viewed as a Temperley-Lieb connectivity. 
This is the last map in Figure 3. In Figure 4, for instance, the perfect matching on the cylinder with 
N = 9 is mapped to a loop configuration corresponding to a connectivity in TL ra= io(0): 



(4.24) 


It is now clear why similarities in the spectra were found in Section 4.1, with the lattice widths of the 
two models related as n = N + 1: The loop segments draw the contours of the web and intersect a 
given slice of the cylinder n = N + 1 times. 

The form of the transfer tangle on the tilted lattice' depends on the parity of n: 



(4.25) 


We call it the tilted transfer tangle. It decomposes as a linear combination of connectivities weighted 
by powers of pi = sin v and p 2 = cosv, where we set 


cos v = 


1 

y/\ + a 2 ’ 


sinu = 


a 

Vl + a 2 ’ 


(4.26) 


corresponding to (4.22) with k = y/l + a 2 . 

Some connectivities contributing to (' T>(v)) M / 2 contain loop segments that will form contractible 
loops when the lattice is folded on the cylinder. The gluing operator (3.13) sends these to zero and 
maps the surviving connectivities to their correct weight, 2 n ' y in this case. This yields 


Z dim = \ (1 + a 2 ) MN/2 F((P(v)) m/2 ) 


(4.27) 


From (3.14), the partition function for dimers can thus be calculated using the transfer matrices of 
critical dense polymers on the tilted lattice, and we arrive at 


n 

Z dim = i (1 + a 2 ) MN/2 ( d + 1 ) Tr p d ((T>(v)) M/2 ). (4.28) 

d= 0 

d=n mod 2 


The rest of the evaluation of Z dim using the tilted transfer tangle is similar to the one of Z cdp on 
the regular (untilted) lattice. One first establishes an inversion identity for T>(v) and then uses it to 
compute the eigenvalues in terms of quantum numbers vy- and rfe. This is done in Section 4.5. As argued 
at the end of that section, the admissibility criterion for a vector (y, r) to appear in the spectrum of 
Pd{T>{y)) turns out to be identical to the similar criterion for pd{D{u , 0)). Partition functions restricted 

6 We note that the factor of | is due to the difference of 1 between n 7 and L. It is ultimately responsible for the relative 
factor of 2 between (2.44) and (3.30). 

7 Note that one can instead use transfer tangles that are reflected vertically compared to (4.25). Then, in Section 4.4, 
arriving at a relation akin to (4.32) would require using the TL„ (0) representation t* defined by r*(ej) = (^(ej)) . 
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to a given d are subsequently obtained by taking the trace of (X>(u )) A/ / 2 in p^, and the steps follow 
those used for the untilted lattice. The conformal, d-specific partition functions are also the same up to 
a redefinition of the nome q. By weighting them by d + 1 and computing their sum, one finally arrives 
at (2.44). 

4.4 Dimer representation of the Temperley-Lieb algebra 

Sections 4.2 and 4.3 established a connection between the dimer model and the Temperley-Lieb algebra 
through the series of maps in Figure 3. This did not involve Lieb’s transfer matrix, but as shown in [23], 
T 2 (a) directly ties the first and last items of Figure 3 via the so-called dimer representation of TL„(0), 

r : TL„(0) End((C 2 )® n_1 ), r(J) = I, r{e 3 ) = aj_,a+ + a+aj +1 , (4.29) 

using the convention cr^ = = 0. We note that T>(v) admits a simple expression in terms of the 

generators of TL n (0): 


Lfl L^J 

T>{v) = (fj \ X 2 j(v) S J, X-j(v) = I cos v + ej sinu. (4.30) 

3 = 1 i =1 

Here onward, we use the variable n (instead of IV + 1) in computations in the dimer model. The 
transfer matrix squared T 2 {a) is, up to a constant, the matrix representative of X>(u) in the dimer 
representation. Indeed, by rearranging the terms of the first exponential in (2.15), we find 

Lfl L^J 

r2 («) = n (l + a (a 2 -._ 2 + <t+_, <r^-)) (j + a ( a 2j -1 a tj + a 2j+ 1)) ( 4 - 31 ) 

3 =1 3 = 1 

as well as 

T ( X ») = \t+a*)V* ' Tmv)) = (1 + a 2 )^-D/2 ’ (a = tanu). (4.32) 


4.5 Lattice integrability and integrals of motion 

Based on examples for small n, it is not hard to see that the commutator \D(u),T>(v)\ is non-zero 
in general. This reflects the non-commutativity of T 2 (a) and T 2 (c/) in the dimer model. As we now 
show, the description of T 2 (a) in terms of TL n (0) reveals the integrability of the dimer model on the 
lattice. 

Following [67,68] , we fix the spectral parameter and inhomogeneities to have the specific values 


„ _ V t — t — (_V V _V V \ 

u — 2 ’ S — Su — V 9 ’ 2 ’ 2 1 2 > ' ' ' 


in which case D(u,£) becomes 


D(IZ V ) = — 

sin v 


V 


V 

r 




J 1 

* 

i_ 

J r 

V 

^ _ 




By using the diagrammatic identity 


& =sm " & 


(4.33) 


(4.34) 


(4.35) 
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and deforming the diagram in (4.34) in such a way that each of the face operators 
oriented as we see that 

D(lt v ) = T>{v). 


V 


and 


V 


is 


(4.36) 


It is now clear from (4.36) and the discussion below (3.17) why T>(u) and T>(y) do not commute: 
They simply correspond to different choices of the inhomogeneities, f u and However, for any fixed 
v E C, we can construct a one-parameter family of transfer tangles, 


D{u,£ v ) 



(«gC), 


which does commute with T>(v), see (3.17): 


(4.37) 


l‘D(v),D(u,£ v )] = [D(lCv),D(u,Cv)]=0. (4.38) 

Using (4.32), it follows that r(Z?( «,£„)) constitutes a one-parameter family of matrices commuting 
with T 2 (a). Tangles that are independent of u and commute with D(^,£ v ) are obtained from the 
coefficients in the expansion of D(u,£ v ) in powers of sin 2 u, 


n— 1 


D(u,£ v ) = ^ J m (sin2u) m , 


(4.39) 


m =0 


and are here called lattice integrals of motion. For n = 3, for example, Jo, Ji and J 2 are readily read 
off from 


D (u,£ v ) = 4(5 + 3 cos 2v) I + sin2u cos v (ei + e^ + 4 shiv (eie 2 — e 2 ei) 


(4.40) 


+ (sin2v) 2 - \l + \(e ie 2 + e 2 ei) 


The matrix realisation of J m = J m {v) in the representation r commutes with T 2 (a) for a = tanu: 

[ T {Jm(v)),T 2 ( tanu)]=0, m = 0,1,..., n — 1. (4.41) 


As indicated, the lattice integrals of motion generally depend on v, or equivalently on a. 

In the end, lattice integrability is built into the dimer model, albeit in a somewhat peculiar way. 
On the finite lattice, different values of a belong to different integrable families, simply because such 
families do not have the same (a-dependent) lattice integrals of motion. In Sections 4.6 and 5.4, we 
will compute the eigenvalues of D(u,£ v ) and those of the lattice integrals of motion and compare them 
with the spectra of the conformal integrals of motion for c = —2. We will find that the a-dependence 
nicely factorises and that the conformal integrals of motion are universal and independent of a, as 
expected. 


4.6 Inversion identities 

From (3.19), the transfer tangle D(u,£ v ) satisfies the following inversion identity: 


£>(«,£„)£>(« + §,&) = I 


( cos(u — |) cos (u + |)) n — (sin(u — |) sin(u + |)) r 
cos (u — |) cos (u + f) — sin(it — |) sin(tx + |) 


(4.42) 
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The function on the right-hand side factorises, yielding the following functional relations for the eigen¬ 
values: 


Eig(D(u,£ v ))Eig(D(u + §,£„)) = < 


n— 1 
2 


2 2n ~ 


kn 


3 = 1 


COS 2 U 


2 / (2j —1 )tt \ 
l 2n / 


sin 


2 • 2 
+ sin v — sin 2u 


2 2 2 
n cos t; 


22n—2 


n 

j=i 


cos 2 u 


-i 2 


sm 


2 (£) 


• 2 *2 
+ sin v — sm 2u 


n odd, 


n even. 


(4.43) 


The general solution is 


L^J 


Eig (D(u,£ v )) = (cosv) n 1 y 1 + a 2 sin 2 1 j + €j sin t 

3 =i 


sin 2 u 


cos v J 


1 + a 2 sin 2 tj + jij sin t 


sin 2u 


cosu J 


(4.44) 

where €j, fij € {+1, —1} and a = tanu, while tj is defined in (3.21). This expression can be rearranged 
to resemble (3.22) and (2.29), 


n —2 


Eig (D{u,£ v )) = (cosv) 71 1 Yl 


, 9 sin 2 2 u \ 9 

1 + ( or -) sm 2 q k 


k =l L 


cos 2 v 


V 1 + « 2 sin2 Qk + 7(^7 sin q k 


y/l 4- /t/ 2 -1-2 sin2 m 


a: sm (/fc 


smft 


l— 


k=n mod 2 

(4.45) 

where = 2 ^ iL , T k = and + r k S {0,1, 2}, while q k is defined in (3.22). 

Remarkably, the transfer tangle D(u, £ v ) provides a unified transfer matrix description of the dimer 
model and critical dense polymers. For a = v = 0, it becomes the transfer tangle D(u, 0) discussed 
in Section 3.6, so (4.45) yields the eigenvalues (3.22) under this specialisation. For u = §, D(u. £„) 
becomes the tilted transfer tangle T>(v ), and (4.45) indeed reduces to the spectra (2.29) of T 2 (a) up 
to an overall factor of = (cosn)”^ 1 , as it should according to (4.32). 

The eigenvalues (4.45) are continuous functions of the parameters u and v, so in any given rep¬ 
resentation, the admissibility criterion for the allowed (z/, r) vectors is independent of u and v. This 
allows us to deduce the following criterion for T)(v) in p^. The constraints on {v, r) are the same as 
those for D(u, 0) given in (3.24). As discussed at the end of Section 4.3, the traces of Pd((^ , (^)) M ^ 2 ) 
in (4.28) can then be computed explicitly, reproducing the expressions (2.44) for the dimer partition 
functions. 

Finally, the analysis of finite-size corrections in Section 4.1 can be extended to D(u,£ v ) simply by 
replacing u(t) in (4.2) and (4.9) by 


u(t) = log (Vl + a 2 sin 2 t + \J\ + a 2 sin 2 u sin t). (4.46) 

Indeed, this function reproduces both w dim (f) and w cdp (f) in (4.3) under the respective specialisations. 
We will return to this in Section 5.4. 


5 Conformal data 

A key objective of this section is to refine the arguments tying the dimer model to a conformal field 
theory with central charge c = —2. Section 5.1 discusses the continuum scaling limit of the dimer 
model, defines the dimer Hamiltonian and identifies the first integral of motion Lq — In Section 5.2, 
we construct expressions for the higher Virasoro modes using their lattice approximate realisations as 
Temperley-Lieb tangles. Section 5.3 then discusses the relation with the fermionic be ghost system 
with c = —2. It culminates with Theorem 1, which asserts that, as modules over the Virasoro algebra, 
the dimer Fock spaces decompose into direct sums of Feigin-Fuchs modules and exhibit reducible yet 
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indecomposable structures. Section 5.4 then examines the conformal integrals of motion for Feigin- 
Fuchs modules at c = — 2 and finds complete consistency with the eigenvalue analysis of the lattice 
integrals of motion. Finally, in Section 5.5, we show that the expression for Lq — ^ found in Section 5.1 
is also compatible with a c = 1 realisation of the Virasoro algebra. Accordingly, we explicitly construct 
Virasoro modes for two non-interacting fermions and find in Theorem 2 that, under their action, the 
dimer Fock spaces are completely reducible Virasoro modules. However, the transfer tangle D(u. £,.) 
is found not to be a generating function for the c = 1 integrals of motion. 

5.1 Scaling limit and dimer Hamiltonian 

In Section 2.3, we used a Jordan-Wigner transformation to express T 2 {a) as a product of the mutually 
commuting operators T q , each written in terms of the lattice fermions Q q and (} q . From (2.28), the 
logarithm of T 2 (a) can be written as 

L"J 

- | log T 2 {a) = arcsinh (a cos ( - Cj ( q + Cl- q C n -q) • (5-1) 

q=i 

For n odd, it has a single ground state, 

|H) =CtiCL3 ---C 2 t Ci t |d), (5.2a) 

2 2 

whereas it has two for n even, namely 

m = cUcU • • • cjclld), \n') = ckU • • • clcJld) = dl«>- (5.2b) 

2 2 2 2 2 

Eigenvectors of T 2 (a) are naturally built from the reference state |d) in the form C91C92 ■ • ■ Cgr|d)- 
To study the scaling limit, however, it is convenient to use a basis for the fermions that builds excitations 
on the ground state(s) instead. We thus proceed by introducing 

= = <n +q (0 <<?<§), (5.3) 

which merely amounts to a relabelling. Here, the mode labels q of Tq and Tq, respectively, take integer 
and half-integer values for n even and odd. The non-trivial anticommutation relations are 

{T p , T*} = 5 p>q , {T p , T^} = d P! g. (5.4) 

From (5.1), we have 

— \\ogT 2 {a) = ^2 arcsinh (a sin ^)($J$ ? + — l), (5.5) 

0 ^9^^ 

where the primed sum indicates that the index q runs over half-integer values for n odd and integer 
ones for n even. For q > 0, the modes Tq and Tq annihilate the ground states (5.2), so the ground-state 
energy can be directly read off: 

Eq = — ^2 arcsinh (a sin ^). (5.6) 

0 <9^^ 

As will be discussed in Section 5.5, T g and Tq can be seen as non-interacting fermions. 

We now turn to the continuum scaling limit , or scaling limit for short, of the dimer model. Roughly 
speaking, this limit retains the states that are the lowest excitations over the ground state and discards 
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those whose energy differences with the ground-state energy do not scale the right way as n —>• oo, see 
(•’•*)• 

This can be defined more precisely. The dimer model actually constitutes a family of models 
characterised by their lattice width and height. The transfer matrices act on the vector spaces (C 2 )® n_1 
whose dimensions increase with n. Because some properties of T 2 (a) are specific to n odd and even, it 
is necessary to study the two parities of n separately. One can then define the sequences of eigenstates 

... <j>t iff ij/t . ..ijrt |Q\ ( 57 ) 

pi P2 Pi 91 92 qj\ n 

of T 2 {a) in (C 2 )® n_1 , made of one state for each n (in steps of 2, so the parity remains fixed), and 
study the behaviour of these sequences of states as n —>• oo. In (5.7), it is understood that i and 
j do not vary with n, but that pk and qe may. Here, focus is on sequences of eigenstates of T 2 (a), 
whereas in Section 5.2 we will discuss sequences similar to (5.7), but based on eigenstates of the dimer 
Hamiltonian H defined below. The operator whose eigenstates are used to define the scaling limit, that 
is, T 2 (a) or H, will be referred to as the reference operator. 

A sequence (5.7) will be said to be conformal if the following two conditions are satisfied: 

(i) The energies E(n ) of the sequence satisfy 

lim n (E(n) — Eo(n)) = k for some n < oo, (5-8) 

where Eq(ji) is the ground-state energy for given n; 

(ii) The limits lim n ^. 0O pfc(?r) and lim, woo qfin) exist for all 1 ^ k ^ i and 1 ^ i ^ j. 

As discussed in Section 4.1, the eigenvalues of conformal sequences have finite-size corrections of the 
form (1.1). Each conformal sequence is expected to converge to a state in a Virasoro module, with n in 
(5.8) equal to 7n?(A — Ao). The difference in conformal dimensions is thus finite, and the corresponding 
conformal state is often referred to as a finite excitation. 

A sequence for which at least one of the above conditions fails is not conformal. In particular, if 
the limit in (5.8) diverges, the energy difference with the ground state is too large and the sequence is 
said to drift off to infinity. Such sequences completely decouple in the conformal field theory. 

In this description, the sequences associated with the states |fl) and on increasing 

values of n are two examples of conformal sequences. In contrast, is not a conformal sequence 

(neither condition is met). We say that labels p t or q.j satisfying condition (ii) lie in the vicinity of zero. 
In the <f\ formulation, q(n) must instead lie in the vicinity of % for the sequence Cg( n )l^) to be conformal, 
meaning that lim n _ > . 00 (g(n) — 1}) must exist. We shall denote by (C 2 )®°° the infinite-dimensional vector 
space generated by the scaling limit of the conformal sequences and refer to it as the dimer Fock space. 

Mimicking the translation from the Schrodinger to the Heisenberg picture, the concept of conformal 
sequences is lifted from the states to the operators: A sequence of operators 0(n) € End((C 2 )® n_1 ) 
is conformal if the action of 0(n) on any conformal sequence is also a conformal sequence or a linear 
combination thereof. Conformal sequences of operators survive the scaling limit. In Section 5.2, we 
shall encounter operators of the form 0±(n) + 02 (n.) where 0±(n) is conformal, whereas O 2 (n) is not, 
as it maps any conformal sequence to a non-conformal one. In this scenario, only 0i(n) survives the 
scaling limit. 

In (5.5), the operators and Tg that are conformal are those for which q is near the lower bound 
of the sum. The large n expansion of logT 2 (a) thus starts out as 

- \ log T 2 (oi) = nf hvL Ik + /bd y + — ( - 777 - + ^2 qi&l&q + V’JV’g)) + o(n~ 2 ) (5.9) 

n 935 0 

where c e g is given in (4.8) and 

^9 = *^9 la =0 ’ ^9 ~ ^ 91 o =0 ( 5 . 10 ) 
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In the scaling limit, the 4 term is expected to become L$ — 

lQ g T2 ( a ) - n /bulk - /bdy) L 0 - ^ qtyqtq + V’W ■ ( 5 ‘ 11 ) 

9^0 

As discussed in the introduction, analysing this 4 term is insufficient to determine the central charge. 

For finite n, log T 2 (a) and log T 2 (c/) do not commute, because the fermions 4> g and in (5.5) 
depend on a , and the relations (5.4) generally hold only if both entries of a given anticommutator are 
evaluated at the same value of a. In the scaling limit, however, the a dependence of the 4 term takes 
the form of an overall factor of a, so in this limit, logT 2 (a) and logT 2 (c/) commute at order 4. 

The dimer Hamiltonian is defined as 

n —2 

=- K + a J+ 1 + a 7 a i+ 1) ( 5 - 12 ) 

j =i 

and can be written in terms of the Hamiltonian tangle H of the Temperley-Lieb algebra as 

n— 1 

H = t(H), « = -^e i GTL n (0). (5.13) 

3 = 1 


H = — -j— log T 2 (a) 


The Hamiltonian tangle is also the leading non-trivial term in the Taylor series expansions of D(u, 0) 
and T>(v): 

D(u,0) = I-2uH + 0(u 2 ), T>(v) = I -vH + 0{v 2 ). 

Expressing H in terms of the fermions is done using (5.5) and (5.12) and yields 


H = 2 Y si n (^)(^g < / , g + ^ q ) + 

Its 4 expansion has the form 


1 — cot(^), neven, 
nodd. 


^ 2 n 

1 — CSc(t^), 


rr_ 2n , i , 271 ( C eS 


+ Y q +^W) +°( n 2 )- 


9^0 


(5.14) 


(5.15) 


(5.16) 


Up to the multiplicative factor of 27r, the 4 coefficient is exactly Lq — given in (5.11) 


5.2 Higher Virasoro modes at c = —2 

One direct way to extract the central charge of a critical model is to construct the Virasoro modes 
explicitly. For lattice models with underlying representations of the Temperley-Lieb algebra, this can 
be achieved via the so-called Virasoro mode approximations. The original idea is due to Koo and 
Saleur [47] and was more recently applied to the gl( 111) spin chain [69,70]. First, the Hamiltonian 
densities are defined as 

n—1 

n k = ~Y cos ( Z T i ) e 3 e T M°)> k € Z^o- (5.17) 

3 = 1 

For k = 0, this is just the Temperley-Lieb Hamiltonian, which in the scaling limit is believed to converge 
to Lq — up to a multiplicative scalar and addition of a constant term. For k > 0, in the scaling limit, 
Hk should converge to the linear combination L *. + L_^, after an appropriate rescaling. In fact, this is 
believed to hold for any TL n (0) representation, in particular in the dimer representation. Because 

[L k + L-ki Lq] =k(L k — L_ k ), (5.18) 
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the commutator [Hk,'Ho\ should then become k (L k — L-k) hr the scaling limit, thus allowing one to 
extract each L k individually. 

Our objective is to investigate this for the dimer representation r, so we apply the Jordan-Wigner 
transformation to 

H k = r(Uk) (5.19) 

and consider the scaling limit. As discussed in Section 5.1, this limit is defined in terms of the energies 
E(n) of a reference operator. In the present context, there are two natural choices for this operator: 
— log T 2 (a) and H = Hq. For finite n, — log T 2 (a) and H do not commute and have different eigenbases. 
For finite n, the expression for H k will of course depend on whether it is written in terms of d> g , 'Fg or 
of <f>q, ipq. Carrying out both computations, however, we find that, in the scaling limit, the two choices 
actually produce the same result at order This is due to the fact that, for q in the vicinity of zero, 

(5.20) 

In the end, the L k Iul are written exclusively in terms of the fermions (j) q and and they are independent 
of the choice of — logT 2 (a) or H as the reference operator. For Lq 1 ™, this is already seen from (5.11) 
and (5.16). We now present the steps of the (slightly shorter) calculation that uses H as the reference 
operator. 

By applying (2.19) and subsequently (2.23), we find 


n— 1 


H k = - 


3 = 1 


n —2 

E 

3 = 1 


Y cos (^) r ( e i) = -Y [ cos ^rr) c } c j+ 1 + cos ( 


n— 1 

llJ ^ 1 )C} +1 C j \ = - Y, Cp, q r)tv q 

p,q=l 


(5.21) 


where 


Cp,q COs( n ) (^p—q— k mod 2n T hp—q,k mod 2n &p+q,—k mod 2n &p+q,k mod 2n) • 


(5.22) 


To take the scaling limit, we restrict k to values much smaller than n. Many contributions to H k 
involve operators that belong to sequences that are not conformal. For instance, for the 5 p+q _k+2n 
contribution, the operator 'i'i p rj q = r]^_ q _ k+2n r]q acts on any conformal sequence by outputting one that 
is not conformal: If q is in the vicinity of p is not. The same applies for the 5p+<j,fc contribution. In 
fact, the last two terms in (5.22) only yield operator sequences that are not conformal. For the same 
reason, the contributions S p - q -k± 2 n and 5 p - q ,k±2n of the first two terms in (5.22) also drop out. The 
only surviving contributions are 


n —1 n —1 

Ak — — ^ ^ VpVq COS ("n‘)^p—<?,—~ ^ ^ VpVq cos {^‘)^p—q,k- 

p= 1 p= 1 


(5.23) 


For k = 0, the results of the computation are already given in (5.15) and (5.16). Setting k > 0, we find 
that, in terms of c^q^q, Ak and Bk are expressed as 


A k = Y' sin (w)^Vfc + Y! M^)h- q %+ Y! sin(^)V>J_ fe </>g, 

O^q^k 


Osiq^^-k 


k<qtZ^ 


B k = Y sin (f )4>\^ q+k + Y' sin (f)^W- Q + Y' sin (f)^_ fc ^> 

0 ^ q < k 


(k > 0). (5.24) 




0 


These expressions for A k and B k can be recast in a simpler form after an additional relabelling, 


</>g q> 0, 
Y q q < 0, 



q > 0, 

q^ o, 


9 


9 


31 


(q = ^ mod 1), 


(5.25) 




where 


(a,/3 € {+, —}). 


(5.26) 


{QpiOq} = fia+/3fl5p+ q) 0, 

Indeed, with this relabelling, we hnd that 

O^q^^-k 

(k > 0). (5.27) 

B k = ^ sin (f) 9 -q 9 t-k + S' sh O^,-A’ 

O^gsS^ OsJg<^-fc 

The scaling limit of H k then reads 

R Hi K‘- + L^=Y! q (ez t e* +t + et q+t e- + ez t e*_ t + e-!: q _ t e-), (*> o). (5.28) 

g^O 

By using (5.11), (5.18) and (5.28), we can isolate each L)( im . In the scaling limit, we hnd that A k 
and B k become L^, im and L \]_™, respectively, with 


rr=E'«(«: AA + A+A). (kez\{o))- (5.29) 

0 

Finally, computing the commutator [L’^ m ,L d ^] = 2kL^ m + -^(k 3 — k), we arrive at 


^o im = E 'q(oz q e+ + e+ q o;)+A, 

g^O 


A = 


0 n even, 
— | n odd. 


(5.30) 


and, crucially, 

c = -2. 


(5.31) 


The weights A = 0 and A = — ^ are readily recognised as the conformal dimensions Agi and A i 2 of 
the c = —2 Kac table, see (4.15). The final compact form for the Virasoro modes is 


Lr = E' 9 (°-X+k + d -q + k0q) + 4,0 A, (k€Z). (5.32) 

g^O 


5.3 Module decompositions 

This section investigates the structure of the dimer Fock spaces as Virasoro modules (under the action 
of (5.32)), comparing them with the representation theory of other realisations of the Virasoro algebra 
at c = —2. 

From (5.32), we hnd that the modes of the helds 0~(z) and 0 + (z) satisfy 

\L% m , 0-} = -{k + q) Ok+gi \ L k m , K] = ~Q 0t +q , (5-33) 

demonstrating that they are primary helds of conformal weight 0 and 1, respectively. In fact, the 
expression (5.32) is a known realisation of the Virasoro algebra at c = —2 and has previously appeared 
in a fermionic be ghost system [71]. It was studied in detail by Kausch [72], 

As discussed in Section 2.4, for finite n, the full state space (C 2 )® n_1 is naturally divided into 
T 2 (a)-invariant subspaces E^_ x labeled by the eigenvalues v of S z . The space E^ l _ l is also invariant 
under the action of the densities 'H k (since it is left invariant by the Temperley-Lieb generators) and is 
thus expected to become a Virasoro module in the scaling limit. In that limit, S z is written in terms 
of the modes as 

S Z = -E' 9 -q 0 t + E 7 °-A + ¥n,0 mod 2 • (5-34) 

g^O g>0 
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Each invariant subspace E^,_ 1 becomes infinite-dimensional, 


K-i EV = span{0 pi • • • 9 TH 9+ ■■■9+ |H)} with j -i = 


v nodd, 
v + \ n even, 


(5.35) 


where the labels pk and (p are non-positive and take half-integer and integer values for n odd and 
even, respectively. A basis of states is obtained by restricting the labels to p\ < P 2 < ■ ■ ■ < Pi and 
Qi < Q 2 < ‘ ‘ ‘ < Qj- The (unique) ground state |fi„) of E v is 


\n v ) 


n' 

0<q<—v 

n' 

0 ^q<v 


V ^ 0, 


v^O, 


(5.36) 


and has conformal dimension 

Here, the indices of the primed products 
respectively. Recalling that 


A„ = 


4u 2 - 1 


(5.37) 


run over half-integers and integers for n odd and even, 


6»+|H)=0 (q > 0) and 0“|fi)=O (q ^ 0), (5.38) 

and comparing with (5.2), we see that the ground state for n odd is |f2) = |flo), while for n even, 
|H) = |fi_ 1/2 ) and \Q'} = |H 1/2 ) = 0 q |fi). 

Because the character over E v is a Verrna module character, the number of states at level l is 
equal to the number p(£) of integer partitions of i. These states, however, are not all descendants of 
|flu). As the following theorem states, E v is not a Verrna module but instead a Feigin-Fuchs module. 
The definition of Feigin-Fuchs modules is reviewed in Appendix A. 


Theorem 1 The module E v is isomorphic to the c = — 2 Feigin-Fuchs module J-\ with A = | — v. 
Thus, its structure is 


(i) for n odd, 


(ii) for n even, 


E v ~ E~ v 


E v : 


A„ 


- O - 


A„ A v+2 A v _|_4 A„ + g 

• ©•©•©• 


A v -|-2 A ^_ j _3 A 1j _(_4 A ^_ j _5 

* • --o-> • ^-o - 


A — t ; A _ t> _{_2 A _^_ j _4 A _^_ j _5 

O -*- • -- O -*- • - O -* • — 


(v > 0), (5.39) 


(v > 3), 

(5.40) 

(v < 


Here, each circle represents an irreducible composition factor whose highest weight is indicated just 
above it. A proof of Theorem 1 is given in Appendix B. Notably, the structure of E v , in the form 
of chains of irreducible modules with alternating arrows, is similar to the structure of Ef l _ 1 as a 
Temperley-Lieb module, the difference being that the chains in the latter case are of finite length [23]. 
We also note that Lq im is diagonalisable on each E v , as are the transfer matrix and Hamiltonian on 
Ef_ l for finite n. 

The structures of the Virasoro modules underlying the scaling limit of the standard modules 
for critical dense polymers are simpler. They contain either one or two composition factors, see (4.16). 
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If there is a single composition factor, that is if n is odd, or if n is even and d = 0, then the module is 
irreducible. For n even and d > 0, there are two composition factors and the scaling limit of was 
conjectured in [73] to be an indecomposable highest-weight module, with extra evidence given in [74]. 
Remarkably, this is actually a consequence of Theorem 1, as shown at the end of Appendix B. 

We conclude this section by recalling the relation between the c = —2 realisation of the fermionic 
be ghost system and that of symplectic fermions [75]. In the chiral case, the symplectic fermions x + ( z ) 
and x~( z ) are two interacting fields defined as 

X a (z) = ]T X?*-*- 1 , (« € {+, -}), (5.41) 

q£l+r 

whose modes satisfy the anticommutation relations 

{Xp,Xq} =pd al3 6 p+q! o, (a, 13 € {+,—}). (5.42) 

Here, d a & is an anti-symmetric tensor which, for simplicity, we set to 

(£ H - (-°! J) ■ <-> 


In the Neveu-Schwarz and Ramond sectors, the value of r in (5.41) is, respectively, 0 and 
corresponding to integer and half-integer labels. The Virasoro modes take the form 


l 

2 ’ 


T SF 
Liu 


X- q Xq +k - X-g+kXq + 4,0 A, 

q ^0 5>0 



0 Neveu-Schwarz, 
— | Ramond, 


(5.44) 


where the primed sums run over integer or half-integer values, depending on the sector. It is noteworthy 
that the generator Lq f admits non-trivial Jordan blocks [75]. 

One can define a map 

Xq^Oq, Xq'-^-qOq, (5-45) 

which, respectively, maps the Neveu-Schwarz and Ramond sectors to the even and odd sectors of the 
be system described in Section 5.2. Indeed, the commutation relations (5.42) become (5.26), and the 
Virasoro modes (5.44) are mapped to (5.32). 

In the Ramond sector, q is never zero, so (5.45) is bijective, and the two theories have the same 
representation content. In the Neveu-Schwarz sector, however, the generator Xq is mapped to zero in 
(5.45), so the map is not bijective. As observed in [72], in this case, the be system is a quotient of the 
symplectic fermion system, and the representation contents are thus different, with the Jordan blocks 
of Lff not transferred to the Lq mode of the be system. 


5.4 Conformal integrals of motion 

The tangle D(u,£ v ) discussed in Section 4.5 is important for two reasons. First, it uniformises the 
treatment of critical dense polymers and that of the dimer model, in the sense that it reproduces 
their corresponding double-row transfer tangles upon specialising the parameters to v = 0 and u = |, 
respectively (although the modules on which the two tangles act are different). Second, and more 
important for the integrability of the dimer model, D(u,£ v ) provides a one-parameter family of self- 
commuting tangles which also commute with the tilted transfer tangle, 

[D(u,£ v ),D(u',£ v )\ = 0, [D(u,€ v ),T>(v)] = 0. (5.46) 

As discussed in Section 4.5, the coefficients in the expansion of D(u,£ v ) in powers of sin 2 u are lattice 
integrals of motion, namely, operators commuting with the transfer tangle and being in involution. At 
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finite size, there is only a finite number of such linearly independent operators. As the simple example 
n = 3 in (4.40) indicates, their explicit computation is complicated already for small n. This is also 
evident from the expressions given in [76]. 

At the level of spectra, on the other hand, a complete picture can be obtained from the eigenvalues 
of the transfer tangle. From Section 4.6, the eigenvalues of —^ log D(u,£ v ) take the following general 
form, valid in any representation, 


n— 2 

Eig( — ^ log D(u, £„)) = — \ log (cos"^ 1 v) + Y ( - uj(q k ) + \(u k + r k ) (u(q k ) - u(-q k ))^ , (5.47) 

k =i 

k=n mod 2 

with the function oj(t) given in (4.46), namely, 

u(t) = log (\/ 1 + a 2 sin 2 1 + \J\ + a 2 sin 2u sin t). (5.48) 

For finite size, the eigenvalues of — 4 log D(u,^ v ) depend on n. Their expansion in inverse powers 
of n has the form (1.1) and the coefficients have been worked out in Section 4.1 for a generic function 
u:(t). We rewrite the expansion as 

OO c~\T) 2)9_1 \ ( \ 

Eig(-^ lo g D(u,£ v )) = n/buik + /bdy + m- )! - hp-i{^ T ), (5-49) 

where the numerical coefficients I- 2 p -1 (^, r) depend on the specific eigenvalue determined by the choice 
of i'fc,Tfc, and, from (4.11), are given by 


h P -i{v, t) 


1 

2 p 


B 2 P (r) + 2 p ^2 {vk + Tfc) 

k^l 

k=n mod 2 


(I)" 1 


(5.50) 


In the present case, the functions \ 2 p -1 of a and sin 2u are the odd-power coefficients in the 
Taylor expansion in t of the function uj(t) given in (5.48). It is not difficult to see that \ 2 p -1 is an odd 
polynomial in sin2u of degree 2p — 1. The first few functions read (with a = tanu) 


Ai(o;, u) 


sin 2u 
cosv ’ 


A 5 (a, u) = 24 


sin 5 2 u 
cos 5 v 


A 3 (a, u) = 2 


20(3 a 2 + 1) 


sin 3 2 u 
cos 3 v 

sin 3 2 u 
cos 3 v 


- (3a 2 + 1 ) 


sin 2 u 
cosv ’ 


+ (45a 4 + 30a 2 + 1) 


sin 2 u 
cos v 


(5.51a) 

(5.51b) 


Setting 


A 2 p -i(a,u) = &-i ( a ) ( sin2 u) 2m 1 , 


m= 1 

we can reorganise the expansion (5.49) as a series in sin2u: 


(5.52) 


Eig(—4 log £>(«,£„)) = ra/b u ik + /bdy + Y ( sin 2n) 2m 1 

m^l 


E £:!’(«) 

p^m 


2 p n 2p 1 I 2p -i(v,t) 

(2 p)\ n 2 p _1 


(5.53) 


By exponentiating this expression and comparing with (4.39), one obtains the 4 expansion of the 
eigenvalues of the lattice integrals of motion J m . This furnishes an infinite triangular linear system 
relating the I 2 p-i to the J m , from which we can reasonably infer that the quantities I 2p -i are themselves 
eigenvalues of certain integrals of motion. As argued below, for c = — 2 , they reproduce the spectra of 
the conformal integrals of motion l 2 P -i- 
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For the specific case of the dimer model, for which one simply sets 2u = v in (5.49) (and for 
the critical dense polymer as well, in which case the relevant tangle is D(u, 0)), we therefore recover 
what appears to be a generic fact for conformally invariant spectra of transfer matrices: The large 
n eigenvalues involve coefficients I 2 p-i which are eigenvalues of integrals of motion. This general 
observation has been first put forward in [28,29], and later checked explicitly in the tricritical Ising 
model [77], for critical dense polymers [76] and in the Ising model [78,79]. 

The explicit values of l 2 p-i given in (5.50) come directly from lattice computations. In the 
scaling limit, the integrals of motion I 2p _1 are universal conformal operators acting in the conformal 
representations occurring in the model. (This makes them much better behaved and, hence, easier to 
deal with, than their finite-size counterparts J m .) Algebraically, the operators I 2p _1 form an infinite¬ 
dimensional Abelian subalgebra of the enveloping algebra of the Virasoro algebra. Although their 
explicit form is not known for all values of p, they can be constructed as charges associated to normal- 
ordered polynomials of degree p in the stress-energy tensor and its derivatives [30,80]. The first few 
read 


Ii 

Is 

I 5 


L 0 — 


c 

24’ 


2 L- n L n + L 


n =1 


c + 2 c(5c + 22) 

- Lq +- 

12 2 880 


bm+n+pfi '■ L m L n Lp : + ^ £ L\-2nL2n 

m,n,p^7j 


-1 


+ £ 

72—1 


72 — 1 

h±V_£_i 


c(3c + 14) (7c + 68) 
290 304 


c + 4 r2 , (c + 2)(3c + 20) r 
~ L ° + -576- L °’ 


(5.54a) 

(5.54b) 


(5.54c) 


where the normal ordering means that the modes are ordered from left to right by increasing values 
of their indices. Because different triplets ( m,n,p ) can produce the same ordered product, some 
L m L n L p will appear with (integer) prefactors in I5. For example, Lq, L_ 2 L\ and L^iLqLi appear 
with coefficients 1,3 and 6, respectively. 

Our main goal here is to verify that D(u,^ v ) is a generating function for the integrals of motion 
for c = — 2 by comparing the eigenvalues of the first few l 2 P -i on the Feigin-Fuchs modules with the 
results (5.50) coming from the lattice. 

From their explicit expressions above, one is easily convinced that the actual diagonalisation of the 
integrals of motion quickly becomes impractical beyond the first few levels. When they are restricted 
to act on a highest-weight state |va), many more eigenvalues have been computed. For instance, the 
eigenvalues of Ii, I3,..., I15 for |va) are given explicitly in [29] for any value of A and c. In the special 
case of c = — 2, the eigenvalues of all integrals of motion on a highest-weight state are 

1 / 1 \ 4u 2 — 1 

l 2 p— i(vaJ = —B 2 p[\v\ + A„ =---, (5.55) 

where / 2 p-i(v) is the eigenvalue of I 2p _i associated with the eigenvector |v). 

The rest of this subsection is devoted to the comparison of lattice and conformal eigenvalues of the 
integrals of motion for the modules E v , whose structures were given in Theorem 1. Our investigation 
covers all I 2p _1 for the zeroth level, as well as I3 and I5 (Ii is trivial) for the first and second levels. 
We will assume c = —2 in the above expressions for I3 and I5, but will discuss the possibility of having 
c = 1 in Section 5.5. 

Let us note that the expansion (5.49) assumes a definite normalisation for the quantities l 2 P - \ [y, t), 
relative to the normalisation of the conformal integrals of motion of which they are claimed to be the 
eigenvalues. Because this normalisation is independent of the vector (v, r), it is fixed by the ground- 
state analysis. The interpretation of the finite-size corrections to D(u,^ v ) as eigenvalues of conformal 
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integrals of motion is subsequently unambiguous. To avoid confusion, from now on, we will denote by 
l 2 p-i the eigenvalues obtained from the action of I 2 p-i in (5.54) on the module E v , and by /^p— 1 the 
quantities (5.50) obtained from the lattice. These can be easily computed for Fock states constructed 
from fermionic excitations of the vacuum |f2): The vacuum itself corresponds to all 1 'ki T k = 0, while 
the action of 9+ q or 9Z q on the vacuum sets, respectively, v 2q = 1 or r 2q = 1, for q > 0. (The creation 
operator 9$ does not modify the quantum numbers.) The values of <5*. = + T/.. are given for the first 

four excited states in Table 1. As before, we set r = § mod 1, taken in {0, 

Level 0 All modules E v , for v integer or half-integer, have a unique highest-weight state |7To) = 
of conformal dimension A„. From the expression given in (5.36), the state |fi„) corresponds to v 2q = 1 
for 0 < q < v (and all other zero) and to r 2q = 1 for 0 < q < —v (and all other zero). In both cases, 
using the general identity 


Y {am + f3) e = (b £+ i (i + §) - B m (§)), (5-56) 

m =0 


we readily obtain 

j-latt („ \ _ ^ 2 p( r ) 
hp-ii^o) — 2 p 


+ 


P 


2 P~ 


t E' i 2r - 1 


0<g<|u| 


%C+ ^ ' l y^\ m + r) 2p- i = B 2p (H + ^) 
2 p 2 p _1 ^ 2P 

m =0 


(5.57) 


which exactly reproduces the result / 2 p-i(vA„) quoted in (5.55). 


Level 1 For a given v, all states in E v are constructed by keeping the same balance between the 
number of 9+ and 9 q as for |7To), but by allowing larger values for p and q. There is a single state in 
E v with Lo-eigenvalue equal to A„ + 1, given by 


7Ti) 




Ci IT cc u<o ’ 

0<g<— v— 1 

9\9~ k |n) u = 0, 

2 2 

»Vj) n' <AI«> ®>°. 

0^q<v— 1 


(5.58) 


where the symbol ~ means an equality up to a possible multiplicative sign. The quantum numbers 5k 
for 17Ti) are given in Table 1. For r/0, the vector 5k = Vk + ifc corresponding to 17Ti) is obtained from 
that of |7To) by changing ^ 2 \v\-i : 1 0 and $ 2 \ v \+i : 0 —> 1 for v / 0. For v = 0, |7Ti) = 1 / 2 C 

contains two excitations, so the difference with |7To) is that <5i : 0 —» 2. We therefore obtain 


C-iO = 


B 2 p(|u| + ±) + V 


2 p ’ 2 p- 1 

also valid for v = 0. Specialising to p = 2, 3 yields 


CO = A* + 6A,„ + rn, #*(*n) = A 3 + f A 2 + f A, + §§, 


latt / 


+ \f r ~' - (M - F 


,3 , 59 a 2 1 15 


where we have used 


B 4 (x) = x 4 — 2x 3 + x 2 — A-, B 6 (x) = x° — 3x 5 + 4T-A- + 


1 


5x 4 x 2 1 


30 


2 42 


(5.59) 


(5.60) 


(5.61) 
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k 

2\v\ - 3 

2\v\ - 1 

2\v\ + 1 

2|u| 4- 3 

Ko) 

1 

1 

0 

0 

ki) 

1 

0 

1 

0 

|tt 2 ) or k^) 

1 

0 

0 

1 

1 ^2) OT M 

0 

1 

1 

0 


Table 1: The quantum numbers 5f : of the first four states in E v for |u| ^ 3/2. For these states, <5^ = 1 
for k < 2|u| — 3 and 5k = 0 for k > 2\v\ + 3. The cases v = 0, ±^, ±1 are special. For v = 0, the pairs 
(< 5 i,< 53 ) corresponding to the first four states are, in order, (0,0), (2,0), (1,1) and (1,1). For |u| = 
the pairs ( 62 , £ 4 ) are (0,0), (1,0), (0,1) and (1,1). For |u| = 1, the triples (<5i, £ 3 , £ 5 ) of the first four 
states are ( 1 , 0 , 0 ), ( 0 , 1 , 0 ), ( 0 , 0 , 1 ) and ( 2 , 1 , 0 ). 


On the conformal side, the module structures described in Theorem 1 show that | 7 Ti) is a Virasoro 
descendant of 1 7To) except if v = as confirmed by a direct computation giving L_i| 7 To) = (v — 5 )|vri) 
and Ti| 7 Ti) = (v + ^) |tto) - I n all cases, we obtain 

L-iLi\iri) = 2A„|7Ti), L_iL 0 £iki) = 2A 2 |7Ti). (5.62) 

The calculation of the eigenvalues of I 3 and I 5 on 1 7 Ti) is straightforward from (5.54) and reproduces 
the expressions (5.60) exactly. 


Level 2 A basis for the two states ^ 2 ), l^) at level 2 in E v is 


M = 9 + _ {v+l) 9 v _, Ivro), |tt'> = 0 v —3 9 + _ {v+h) |vr 0 ). (5.63) 

It follows that for v ^ 0, ^ 2 ) differs from | 7 To) by one new excitation with mode label k = |u| + ^ and 
none for k = |w| — §, whereas 1^) has an excitation for k = H + § and none for k = \v\ — For 
v ^ 0 , the two mode labels for 17 r 2 ) and 1 ^) are identical but interchanged between the creation and 
annihilation operators. Thus, 


For p 
and 


2,3, 


4p- 1(^2 °r v t' 2 ) = 

4p -M or vr 2 ) = 
and in terms of A„, 
4 att (7T2 or ir' 2 ) = 

or 7 r') = 


B2g(l 2 / + a + £t((H +IW* - (kl - I) 2 " -1 ). 
B2,,( 2 l + " ) + £t((H + if’- 1 ~ (kl - I) 2 ”"). 

we obtain the following two pairs, 

A 2 + 12A V + ff ± 3^8A V + 1 

A v + —+ 60A^ + ± 15(A tl + ^)\/ 8 A^ + 1. 


(5.64a) 

(5.64b) 

(5.65) 

(5.66) 


Let us now compare these values with the eigenvalues of I 3 and I 5 on the two states at level 2 in the 
modules E v . According to Theorem 1, these two states are descendants of |7To) except for v = 0, |, §. 
It is not difficult to compute the action of the two integrals of motion on the basis L_ 2 | 7 To), L 2 _ 1 | 7 To). 
We find (still with c = — 2) 


level 2 


A 2 4- 12A 4- 


+ 


-2 12A, 
6 4 


level 2 


A 3 + 

*-*1) I A 


+ 


( 45A„ + 60A 2 + 45A„ 

V 30A 1; + f 75A„ + ^ 


(5.67a) 

(5.67b) 
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In each case, a simple diagonalisation shows that the eigenvalues are precisely the values given in (5.65) 
and (5.66). 

For v = 0, the highest-weight state 17To) has weight Ao = —|. In this case, the eigenvalues in (5.64) 
are degenerate (for all p ). From the module structure in (5.39), the state |7To) has a single descendant 
at level 2 , as follows from the identity L_ 2 \tto) = 2 L 2 _ 1 | 7 To). This single descendant is an eigenstate of 
13 and I 5 and accounts for one of the two degenerate eigenvalues, namely 

T — A 2 _l19A i 599 _ 3 367 t _ a 3 , 119 A 2 , Rn A i 4 285 _ 61457 /cr R o\ 

-*3 — Aq + I2Ao + 12 0 — 960 ’ -*5 — ^0 T 4 Ag + b(JAo + 336 — 4Q752- (5.6oJ 

The other state at level 2 is not a descendant of 1 7To) . Instead, it is a highest-weight state with weight 
A 2 = ir. As such, it must be an eigenstate of all l 2 P -i with eigenvalues given by (5.55). For p = 2, 3, 
these eigenvalues are 


/ 3 (va 2 ) = 1B 4 (|) = 


3 367 
960 ’ 


4 (va 2 ) = |b 6 (§) = 


61457 
10 752’ 


(5.69) 


as expected. 

For v = 5 , A 4 / 2 = 0 and the relevant quantities, (5.65) and (5.66), are equal to 


r datt 5 flatten r 239 959l lrlatt / '„ \ rlatt^/n r 505 8065 1 

ih ( 7r 2),I 3 \^2)t - il20’ 120 )’ i J 5 {n 2 ),l 5 (7T 2 )> - joog, -oQg-j. 


(5.70) 


Because L_i| 7 To) = 0, the state 17To) has a unique descendant at level 2, namely L_ 2 | 7 To). The action of 
I 3 and I 5 on L_ 2 ( 710 ) is most easily obtained from the upper-left entry in the 2x2 matrices given in 
(5.67a) and (5.67b), namely 


I 3 — + 12A„ + 


_ 239 


a 4 


The second state at level 2 is the descendant L_i| 7 Ti) of the highest-weight state 1 7 Ti), which was 
encountered at level 1 and has A 3 / 2 = 1. The eigenvalues of I 3 and I5 on such a state are given in 
(5.60). For A 3 / 2 = 1, they become 


h — + 6 A t , + — 


959 

120 ’ 


h = K + + 


15 A 

2 1 


I 253 _ 
_r 336 


8 065 
336 ’ 


(5.72) 


once again in full agreement. 

For the last exceptional case, v = |, (5.65) and (5.66) give 




3 239 1079 1 
120 ’ 120 J’ 




r 61 237 8 317 i 
1 336 ’ 336 J 


(5.73) 


A direct computation, using the expressions (5.54) and the basis {#ji 3 0g"|f}), shows that 

the first is an eigenvector of both I3 and I5, and that the eigenvalues indeed coincide with (5.73). 

We thus find that the c = — 2 modules in Theorem 1 are entirely consistent with the first finite- 
size correction terms of the eigenvalues of the dimer transfer matrix (squared), related through the 
eigenvalues of universal integrals of motion acting on these modules. We note that the finite-size 
corrections for critical dense polymers are given by the same numbers I 2 p -1 at all orders, albeit generally 
with different degeneracies. The Virasoro modules in critical dense polymers (discussed in Section 5.3 
and Appendix B) are structurally different [73,74] from the ones in the dimer model. We note, however, 
that the eigenvalues of the integrals of motion only depend on the composition factors present and are 
independent of the embedding structures. The consistency with c = —2 can therefore also be confirmed 
in this case [76]. 


5.5 Virasoro modes at c = 1 

In (5.11), we gave an expression for Lg im — ^ in terms of the fermions <f> q and ip q and built a c = —2 
realisation of the other Virasoro modes in Section 5.2. In this section, we show that Tg im is also 
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consistent with a c = 1 conformal description, and analyse the resulting Virasoro modules. However, 
we also argue that this c = 1 picture is irreconcilable with the eigenvalues of the integrals of motion 
I l 2 p-\ discussed in Section 5.4. 

Let us start by recalling the fermionic construction of the Virasoro modes for c = A single 
fermionic field has a mode expansion of the form 

f{z) = E with {fpJq} = S P+q,0 , (5.74) 

q 


where the modes have integer or half-integer labels, depending on whether we consider the Rarnond or 
Neveu-Schwarz sector, respectively (the primed sums run over the corresponding values). The Virasoro 
modes take the form 


L k = Y^ti-^fk-qfq + X 1 4 , 0 , 
q> I 


A lf = 


^ Rarnond, 

0 Neveu — Schwarz, 


(5.75) 


and satisfy the Virasoro commutation rules for c = By computing [L^,f q ], one confirms that f(z) 
is a primary field of conformal dimension A = 7. 

The fermions used in the previous sections satisfy {#+, d~} = S p+q and give rise to two independent 
fermions of the type given in (5.74), defined by 



6 t 

Q> 0, 


0 q q > 0, 



II 

o 

to 

II 

7 #(l^o) 9 = 0, 

(5.76) 

Vq 

q< o, 


0+ q < 0, 



in terms of which the anticommutation rules read 


{ft\f^} = S a ,pS p+q ,o, («,/?€ {1,2}). 


(5.77) 


The Rarnond and Neveu-Schwarz sectors thus correspond to the even and odd parities of n of the dimer 
model. As the fermions f q 1] and f q 2) are independent, {/p X) , /) 2) } = 0, the sum of their Virasoro modes 
yields a realisation of the Virasoro algebra with c = 1 in the form of two copies of (5.75), 


r|'= 

?>■§ 


+ fklf?) + & K, c 


A 2f = 


| n even, 
0 n odd. 


(5.78) 


In addition, the following operator, 

* = i( E" f-qfq 2 ' - E' f-lfq') ( 5 ‘ 79 ) 

q^O q> 0 

commutes with all the Virasoro modes: 

[Ltx] = o, (fee Z). (5.80) 

Rewriting the modes L^ in terms of the fermions 9± provides the announced conformal realisation 

with c = 1. We also note that for k = 0, Lq — ^ = Lq — ^ given in (5.78) for c = 1 exactly reproduces, 
for the two parities of n, the expression of L Q im — ^ = Lg im + ^ found in (5.30) for the dimer model and 
c = —2. The other modes, for k ^ 0, have different expressions in the two pictures and are expected 
to lead to Virasoro modules with different structures. Crucially, the modes L| f were written down 
heuristically from the observation that Lq 1 " 1 appears to be described by two fermions, and not from the 
Virasoro mode approximations as in Section 5.2. 
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Likewise, the operator X differs from the operator S z used before and therefore partitions the 
dimer Fock space (C 2 )®°° in a different way. In particular, the states |fl„) defined in (5.36) are not 
eigenstates of X for |w| ^ 1. To find an eigenbasis of X, we define two new fermions by setting 

fW = !=(§++§-), /< a > = -±=(0+-0-), (all?) (5.81) 

so that the two sets of modes 9^, 9± are linearly related. One can easily verify that the anticommutators 
are preserved, 

{9p,9^} = 5 a+ /3fiSp +q fi, (a, /3 €{+,—}). (5.82) 

In terms of 9±, the operator X is given by 

x = — ^2 O-qOq + mod 2 (5.83) 

< 3^0 q> 0 

and has exactly the same expression as the operator S z had in terms of 9±. Also, the Virasoro modes 
acquire a rather simple form, 

£”= EL-l)fcA + + »LA') + A ” { M- < 5 - 84 ) 

<?>§ 


Because of 

[I%,e±] = -$+q)0± +q , (5.85) 

the fields 0 ± ( z ) are primary, both of dimension A = ^. 

The states 

8pi---8pA---°i\ n ) ( 5 - 86 ) 

form a basis of the dimer Fock space (C 2 )®°°, which is the direct sum of the subspaces E x where X 
acts as x ■ I. These subspaces are generated by the states of the form (5.86) with j — i = x for n odd 
and j — i = x + ^ for n even. Each E x is invariant under the action of the modes L| f and is thus a 
module over the Virasoro algebra at c = 1. The (unique) ground state |fl x ) of E x is 


\& x ) 


9_ k |fi) x V 0, 

0 <k<—x 

< 

[J' 0+Jfi) xfiO, 
0^k<x 


(5.87) 


and has conformal dimension 



(5.88) 


Comparing with Section 5.2, we find that 9± and 9± realise identical anticommutator algebras. 
However, the Virasoro modes (5.32) and (5.84) have different forms. The next theorem is the counter¬ 
part of Theorem 1 and gives the structure of E x as a c = 1 Virasoro module. 


Theorem 2 The module E x is an irreducible c = 1 Virasoro module, except for n odd and x = 0, in 
which case it is an infinite direct sum of irreducible modules. Thus, its structure is 


E 


X 


' A 2 0 f A 2 2 f A 2 4 f A 2 6 f 

• ©•©•©• 

A 2 , f 


n odd, x = 0, 


otherwise. 


(5.89) 
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A proof of the theorem is presented in Appendix C. 

As already indicated, a comparison of the lattice and conformal integrals of motion involves relative 
normalisations. For the values cq € { — 2,1} of the central charge c, the question is thus whether the 
lattice results m atch the spectra of ^p-i 1 ' 2 p-i | c=co f° r some normalisation coefficients ^p-v I n 

Section 5.4, we presented evidence that the c = —2 integrals of motion reproduce the lattice results I 2 p -1 
if ^ 2 p-i = 1 f° r P- Thus, to match f° r given p, the corresponding spectrum of 7 ^, l 2 p-i | c=1 

following from Theorem 2 must match the similar spectrum of l 2 p-i | r= _ 2 following from Theorem 1. 
Since 

x 2 /2 v 2 /2 

ch \E x }(q) = S— and ch [E v }(q) = ?—— (5.90) 

7 ( 9 ) p{q) 

are equal for x = v, the first integrals of motion Ii| and Ii| __ 2 have matching spectra (meaning 
that 7 f =1 = 1 ), in accordance with the observation that the effective central charge is the same in the 
two cases. 

However, the spectrum of I 3 1 C=1 does not match the lattice results. To appreciate this, for n odd, 
it suffices to compare the 13 -eigenvalues on the Ii-eigenspaces of eigenvalues — ^ and For both 
central charges, these eigenspaces are spanned by highest-weight vectors. For c = 1, the corresponding 
conformal weights are 0 and respectively, while for c = —2, the weights are — | and |, respectively. 
Having matching spectra thus requires that 


73 1 T 3 (v 0 )| c=1 = T 3 (v_i / 8 )| c= _ 2 , 7^ 1 T 3 (vi / 2 )| c=1 = T 3 (v 3/8 )| c= _ 2 , 


(5.91) 


that is, 


_ 3_^ c=1 

320 '3 
,,c=l 


7 

960’ 


AL ~c=i 
320 7 3 


127 

960’ 


(5.92) 


to which there is no solution for 73 i . For n even, it suffices to consider the spectrum of I 3 on the 


four-dimensional Ii-eigenspace of eigenvalue . For both central charges, I 3 is diagonalisable on that 
subspace. For c = —2, all four eigenvalues are j^, whereas for c = 1, there are two distinct eigenvalues, 
namely and each appearing twice. 

These results for c = 1 are thus incompatible with the corresponding lattice result. We conclude 
that D{u,^ v ), of which T 2 (a) is a specific matrix realisation, is a generating function for the integrals 
of motion for c = — 2 , but not for c = 1 . 


6 Discussion 

Are dimers on the square lattice described by a c = 1 or a c = —2 conformal field theory? Previously, 
the evidence supporting the c = — 2 description was unsatisfactory. This paper is an attempt to 
remedy this. First, we reinforced the connection between the dimer model and critical dense polymers, 
by mapping the former to the latter through a series of maps. The dimer partition function was then 
computed using the loop model and its double-row transfer tangles. Second, we used the description of 
the transfer matrix squared in terms of the Temperley-Lieb algebra TL ra (0) and its dimer representations 
to understand the lattice integrability of the dimer model and construct expressions for all the Virasoro 
modes in the scaling limit. The resulting expressions were recognised as producing a c = — 2 realisation 
of the Virasoro algebra previously encountered in fermionic be ghost systems. As modules over the 
Virasoro algebra, the dimer Fock spaces were found to be of Feigin-Fuchs type. For one parity of the 
lattice width n, the corresponding modules are reducible yet indecomposable, while for the other parity, 
they are completely reducible modules. Third, we found that the lattice integrals of motion generated 
by the transfer matrices have eigenvalues that precisely match those of the conformal integrals of motion 
at c = —2. Altogether, we believe that this constitutes compelling evidence for the soundness of the 
c = —2 description of the dimer model. 

We also investigated a possible c = 1 description of the dimer model, compatible with the expres¬ 
sion for Lq 1111 — obtained from the dimer transfer matrix. These c = 1 modes take the form of the 
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c = -2 

c = 1 

Primary fields 

0 + ,0~ 

0 + ,0~ 

Conserved lattice operator 

S z 

*("> 

Virasoro realisation 

fermionic be ghosts 
(non-unitary) 

two non-interacting fermions 
(unitary) 

Module structures 

completely reducible (n odd) 
reducible yet indecomposable (n even) 

completely reducible 


Table 2: Comparison of the c = —2 and c = 1 descriptions. 


sum of two copies of c = | modes, and their action on the dimer Fock space turns out to be completely 
reducible. This description is unitary. However, the corresponding conformal integrals of motion are 
found not to reproduce the eigenvalues of their lattice counterparts, as the c = —2 ones did. 

From the point of view of Lieb’s transfer matrix approach, the c = —2 description therefore seems 
to be favoured. This conclusion is built on the premise that the ^ expansion of the logarithm of the 
transfer tangles/matrices precisely yields the conformal integrals of motion, as has been observed to 
be true for other models [76-79]. Of course, allowing combinations of the I- 2 p -1 still gives commuting 
elements of the universal enveloping algebra of the Virasoro algebra, and it is not clear at this point if 
this could cure the c = 1 interpretation of the integrals of motion. 

Table 2 draws a comparison between the two conformal descriptions. Their only common feature 
is Lq 1 ™ — T^j-, which, perhaps surprisingly, can be expressed in terms of two different sets of primary 
fields. We are not aware of other physical models for which such a duality exists. It thus seems that 
the central charge depends on the questions one hopes to answer. As discussed earlier, Lieb’s transfer 
matrix approach provides a natural framework to study properties of spanning trees, consistent with 
the c = —2 description. On the other hand, the c = 1 description has been found to be suitable for 
questions involving monomers or heights. To investigate this further, a sensible next step is to search 
for a transfer matrix formulation of the dimer model, which leads to a c = 1 description, and investigate 
the ^ expansion of its eigenvalues. Also, we find that the operator 

L^J 

X(n) (vWn-q + VqVn-q) ~ i 6 n ,0 mod 2 (??t 7 ?n - 5) € End((C 2 ) 0n_1 ) ( 6 . 1 ) 

<7=1 

commutes with T 2 (ot) at finite size and becomes X in the scaling limit. As opposed to S z , A (n) appears 
not to have a simple interpretation when expressed in terms of Cj , C]. or so it again seems 

J • J J 

that the c = —2 description is more natural for Lieb’s transfer matrix approach. Instead, X w may 
have a more natural interpretation in a c = 1 transfer matrix formalism. 

The results of this paper raise several new questions for future work. In addition to the c = 1 
issues above, one may ask whether the Temperley-Lieb connection and the conformal description of the 
dimer model extend to other lattices and/or boundary conditions. For instance, it is not clear how these 
extend if we allow monomers on the boundary, define the model on the hexagonal or the triangular 
lattice, or fold the lattice on a torus. The torus problem addresses the bulk part of the system and 
is expected to involve non-chiral representations. An understanding of this will yield valuable new 
insight into the dimer problem and its description as a conformal field theory. We also note that the 
parity of the cylinder circumference M is necessarily even for the series of maps from perfect matchings 
to polymer configurations to be well defined. A detailed analysis of the integrability and conformal 
properties of the dimer model for M odd remains an open problem. 
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A Feigin-Fuchs modules 

The Feigin-Fuchs modules are a special class of modules over the Virasoro algebra. Their structure was 
first described by Feigin and Fuchs [48], see also the book by Iohara and Koga [81] and the summary 
by Ridout and Wood [82, Appendix A], 

One realisation of the Feigin-Fuchs modules is built on the Heisenberg algebra H. This algebra is 
spanned by modes ak,k E Z, and a central element 1 satisfying the relations 


<r] — k Sk+ifl 1, [afc, 1] — 0. 


(A.l) 


A one-parameter family of highest-weight representations J-\ over H, with A € C, is constructed from 
the action of the universal enveloping algebra of H on the highest-weight state | A), with the rules 


afe|A) =0 (k> 0), 

A basis of states is given by 


a 0 |A) = A|A), 11 A> = | A). 


''' °-fe I A) ; r € Z^o, h> k 2 > ■■■ > k r >0, m, ra 2 , - - -, n r E Z >0 }. 


(A.2) 


(A.3) 


It is not hard to see that J-\ is irreducible as a module over H. Indeed, for each element |w) of the 
basis (A.3), one can construct the word 


TYl\u — CL i 


' a kl a kl 


such that 




|w') = k w <5w )W / |A), = JJn*! (A.4) 


i= 1 


for any |w') of the same basis. One can then get from any state |v) = ^ W a w |w) to an Y other state |v') 
by (i) selecting a |w) for which a w ^ 0, (ii) applying on |v) to arrive at |A), and (iii) descending 

from | A) to |v'). 

The universal enveloping algebra of H carries a one-parameter realisation of the Virasoro algebra, 

11 1 oo 

-hfc = "y ^ a q a k—q — (A: + 1 ) Q CLk (k 7^ 0), Z/q = ^ a 0 T y ^ CL—qCLq Q E C, (A. 5 ) 

<jGZ q =1 

for which the central charge is c = 1 — 3 Q 2 . As a module over the Virasoro algebra, the Fock space J-\ 
is referred to as a Feigin-Fuchs module. Its highest-weight state |A) satisfies Tq F |A) = A^|A) with 


Aa = ^A(A-Q). 


The character of Tx is the same as that of the Verma module Va a of highest weight A^, 

A_c_ 

n x 24 

ch[J\](g) = 


nr=i(wr 


(A.6) 


(A.7) 
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This implies that J-\ and Va a have identical composition factors, but not necessarily the same structure. 
Because A a = A q-\, this also applies to However, chfJ^Kg) = ch[J r Q_A](<?) does not imply 

that J-\ and Tq- a are isomorphic. Setting 


Q = \lzzj(p' -p). A = X r,s = -\l%- ( r ~ !) + \I^ZI ( s ~ !) 


pp 

yields the central charge and conformal dimensions 

, {p' - p ) 2 


2 P 


2 pf 


c = 1-6 


j > 


A = A r , s = 


(p'r — ps) 1 — (p — p') 


/A2 


For 


pp ' 4 pp' 

P-/P € Z >0 , p ^ gcd(p,p') = 1, 


(A. 8 ) 


(A.9) 

(A.10) 


the Feigin-Fuchs module J-\ is irreducible unless Aa is a conformal dimension in the Kac table. The 
special values correspond to A = X TtS (and thus Aa = A rjiS ) with r and s integers, and in this case 
7>,. s = J -a is a reducible Virasoro module. Because A r _|_ PjS _)_p' = A rjS and, hence, J^r+p.s+p' — we 
can restrict r and s to be positive integers. 

The module structure of depends on whether (r, s) is (i) a corner entry of the Kac table, 
meaning that r is a multiple of p and s is a multiple of p', (ii) a boundary entry of the Kac table, 
meaning that either r is a multiple of p or s is a multiple of p' , but not both, or (iii) an interior entry 
of the Kac table, meaning that r is not a multiple of p and s is not a multiple of p'. If p = 1, there are 
no interior entries, and if p = p' = 1 , there are no boundary entries either. 

The different module structures are displayed in Figure 5. There, the composition factors in the 
socle — the maximal completely reducible submodule — are represented by black circles. Quotienting 
by the socle gives a new module. The second socle layer is the socle of this quotient module. Its 
composition factors are represented by gray circles in Figure 5. Repeating the construction defines 
the third socle layer , whose composition factors are represented by white circles. Feigin-Fuchs modules 
never have more than three socle layers. The arrows in Figure 5 describe the action of the Virasoro 
algebra. If an arrow points, say, from a composition factor O to a factor •, vectors in • can be 
reached from vectors in O by the action of the Virasoro algebra, whereas no vector in O can be 
reached from •. 

For corner entries, J> )S is completely reducible and is the direct sum of the composition factors of 
VA r , s - For the other entries of the Kac table, is indecomposable. For boundary entries, F r>s takes 
the form of a chain of composition factors that alternate between the socle and the second socle layer. 
As drawn in Figure 5, there are two possibilities, and corresponds to the one where the composition 
factor of conformal weight A rjS + rs is in the second socle layer. For interior entries, the structure is 
that of a braid, with some arrows pointing in the opposite direction compared to the similar arrows 
in the corresponding braided Verma module. The composition factor with weight A r . jS + rs belongs to 
the third socle layer. Only the top three cases of Figure 5 will come up in Appendices B and C. 


B Module structures at c = —2 

The main objective of this appendix is to give a proof of Theorem 1. Our proof builds a family of 
isomorphisms between the Virasoro dimer modules E v and Feigin-Fuchs modules J-\ at c = —2, for 
which 

(p,p , ) = ( 1,2), Q = 1, A = A r , s =- Ar :S = - J -. (B.l) 

The Kac table for c = —2 consists of corner entries (for s even) and boundary entries (for s odd). The 
module structure of J- TyS is discussed in Appendix A. For in particular, for s even, this structure is 

Ai jS Ai )S+ 4 Ai jS+ 8 Ags+12 

J~i,s — d~\ ,4_s : 0 0 0 0 0 0 0 0 ... ( s ^ 2, even), (B.2) 
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Not in the Kac table : 




Corner entry : 


• ©•©•©•© 


Boundary entry : 


< 

O 


O 


• *-o 

or 

O -* • 


• ^- O 

O -* • 


Interior entry : 


O -* O 


O 


O -- O 


O 


< X X X X X X 


Figure 5: Module structures of Feigin-Fuchs modules. The conformal dimensions of the composition 
factors of J-\ are the same as those of the corresponding Verma module Va a , and grow towards the 
right. 


whereas for s odd, 


F\,s ■ < 


' X,4-s 

Al,6-s 

Al,8-s 

Xio-s 

^1,12-8 

^1,14-s 

w < 


* w * 


> W K 


^l,s 

^l,s+2 

X,s+4 

X,s+6 

Al,s+8 

^l,s+10 

o — 


— o — 


— o — 

—► • «- 


(s ^ 1 , odd), 


(s ^ 3, odd). 


(B.3) 


We note that the periodicity property J- r - )S — < 7 >+i jS +2 yields Xs — J 7 ( 4 - s )/ 2,2 for s even and Xs ~ 

X-*)/2,1 for s odd - 

We now relate these modules to the ones found in Section 5.2 from the dimer model. The realisation 
(5.32) of the Virasoro modes depends on the parity of n. We thus define the (parity-dependent) map 


0 : H —» End((C 2 )®°°) 


(B.4) 


whose action on the generators is 

0 ( 1 ) = I, g(a k ) = J2' d ~ q 0 t+k - T! 9 - q +k 9 q + ¥k,0 Sn,l mod 2- (B.5) 

0 q> 0 

Here, we recall that primed sums indicate that q runs over integers or half-integers for n even or odd, 
respectively. Comparing with (5.34), we see that 

g(a 0 ) = ±-S z . ( B -6) 

One also verifies the commutation relation 

[d(ak),s(ae)\ = k5k+e, o, (B.7) 

thus proving that 0 is an H -homomorphism. It is then straightforward to show that 

0 (afc)|H„) = 0, (k > 0), 0 (a o )|fi„) = ( d - v)\Q v ), (B. 8 ) 

for all v. As a consequence, the map 

(1 : -Fa ->• E v , A = i - v, (B.9) 
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defined as 


(BIO) 


f)(« ■ ■ ■ < W) = B(o h )"‘8(o fa )" 2 • • • SKHS!,,), 
is a homomorphism of H-modules. 

Remarkably, the Virasoro and Heisenberg modes satisfy the commutation relations 

[ L fc F ! 0 <?] = ~qa k+q - ^^4+9,0 1, [ L fc im >0K)] = -9 0(«fc+g) - ^TT^4+9,0- ( B - n ) 

The similarity between these commutators is crucial for the proof of the next two propositions. 


Proposition 1 The map h is bijective. 

Proof Both J-\ and E v are infinite-dimensional vector spaces, graded by the eigenvalues of L g F and 
Lq™, respectively. For E\, a state of the basis (A.3) has grade YH=i ki n i- At grade £, the number of 
states is p(£), the number of integer partitions of £. For E v , the number of states at grade £ is also 
p(C), as is obvious from its character (2.43). 

The grade of a basis state |w) of E\ and that of f)(l w )) £ E v can be calculated using only the 
commutators (B.ll) at k = 0. Because the right-hand sides have the same form, the grades of |w) 
and f)(|w)) are equal, so f] is grade-preserving and thus maps each subspace of E\ of grade l to the 
corresponding subspace of E v , also of grade £. Both subspaces are finite-dimensional, and they have 
equal dimensions. We now show that f) is bijective at each grade. 

At a fixed grade £, the domain and target spaces have equal dimensions, so it suffices to show that 
(l is injective. Suppose there exists a state |v) = X^ w a w|w) in F\ such that f)(l v )) = 0, with the sum 
on w running over the elements of the basis (A.3) with JV =1 kj.rii = i. Then, from (A.4), 

Kw 1 0 (m w )fl(l v )) = K w lf l( m w|v)) = a w |^) = 0 -)• a w = 0, (B.12) 

so the only possibility is that |v) = 0 . ■ 

It follows that 1) is, in fact, an isomorphism of H-modules. The next proposition shows that this 
isomorphism property also holds at the level of the Virasoro modules. 


Proposition 2 The map f] is an isomorphism of Virasoro modules. 

Proof By extending q to the universal enveloping algebra of H in the obvious way, we see that the 
operators g(L^ F ) provide a realisation of the Virasoro algebra with c = —2 in the Fock space (C 2 )®°°. 
We also note, from (B.ll), that 0(Tj( F ) all commute with g(ao) and therefore leave the subspaces E v 
invariant. In fact, 

{g{LF)-Lf*)E v = 0 (B.13) 

holds for all v and all k. The proof of (B.13) uses, as a basis for E v , the elements f)(l w )) with |w) 
of the form (A.3), and is inductive on the grade, with starting point |H„) at grade zero. For k > 0, 
fl(L£ F )|n„) = 0 and Lf m \Q v ) = 0. For k = 0, 


g(Ll F )\n v ) = L£ im \n v ) = A v \n v ) (B.i4) 

with A„ defined in (5.37). For k < 0, we have 

$(L F _ F k )\n v ) = £a w J>(|w», £ill^>=EA"f)(|w>), (B.15) 

w w 

for some a w ,/3 w € C. These coefficients can be calculated using 

a w |H„) = k~ 1 lg(m w ), 0(^)1 |O v ), /3 w \tt v ) = re" 1 [g(m w ), L d Ef ] (B.16) 
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Because [ 0 (I/ FI *.),fl(a g )] = [L^,g(a g )], we also have that [ 0 (T FF fc ), 0 (m w )] = [L^£,g(m w )] for any word 
m w , and therefore, a w = /3 W , completing the proof at grade i = 0 . 

For the higher grades, any basis state |s) of E v in the basis (B.10) takes the form of g(aq)|s / ) for 
some q < 0, implying that the respective grades i and i' of |s) and [s') satisfy l > (!. Then, 


( 9 ( L k F )- L k m ) |s) = ( 0 (Lf)-Lr)0(a g )|s') = [(0(Lf)-L^ m ), 0 (a g )]|s')+g(a g )(0(L| F )-L^ m )|s') =0 

(BT7) 

where the inductive assumption was used at the last equality to set the second term to zero. ■ 

It follows from this proposition that E v ~ E\ for A = | — v. Comparing with (B.l), we find 
that the integer and half-integer values of v correspond to A = Ai iS with s = 2(1 — v) odd and even, 
respectively. The module structure of E v can therefore be directly read off from (B.2) and (B.3), thus 
completing the proof of Theorem 1. 

The structure results for the dimer Fock spaces in the c = — 2 description have interesting reper¬ 
cussions for critical dense polymers. The structures of the Virasoro modules of the scaling limit of 
the modules V[( were conjectured in [73] for logarithmic minimal models CM (1, j/), and more gen¬ 
erally in [74] for CM(p,p'). Critical dense polymers corresponds to CM{ 1,2), and in this case from 
Theorem 1, the limiting structures of can be proved explicitly to be given by 


\/d n^-oo 
v n ' 



n odd, n even with d = 0 , 
n even with d > 0 . 


(B.18) 


To prove this, we first recall that the limiting character of is written in terms of irreducible 
ones in (4.16). For n odd, as well as for n even with d = 0, the limiting structure of V([ is completely 
determined by its character and is irreducible. 

For n even and d > 0, the character is the sum of two irreducible ones, so the Virasoro module has 
two composition factors, and its structure is not completely fixed by the character. Recalling that /3 is 
the fugacity of loops in TL n (/3), an injective intertwiner f n/ i was introduced in [83], for all j3 E R, and 
shown to intertwine pd(H) and a spin-chain Hamiltonian HI £ End((C 2 )® n ^ 1 ). At /3 = 0 , H becomes 
the dimer Hamiltonian H, and one can show that f n ^ actually acts as an intertwiner for each of the 
TL n (0) generators: 


fn,d ' Pdifij) = r v (ej) ■ f n ,d, (v = j = 1,. .., n - 1), (B.19) 

where t v denotes the restricted action of r on the subspace E^_ 1 . 

The injectivity of f n% d was proven in [83] and implies that, as a TL n (0)-module, is a submodule 
of E'n_ x . v = ^=5. This is true for all n and all d = n mod 2. Because the Virasoro mode approximations 
are constructed from TL„(0) tangles, the same inclusion, V n ^ C E^_ l for v = also applies to the 
corresponding Virasoro modules in the scaling limit. Comparing the two characters, we find that the 
structure of E v for n even can be read off directly from (5.40) and is indeed highest weight, as in (B.18). 

C Module structures at c = 1 

This appendix provides a proof of Theorem 2, building on the known structures of Feigin-Fuchs modules 
at c = 1, see Appendix A. The corresponding data is 

M) = ( 1,1), Q = 0, A = A r , s = ^, A r , s = ^^. (C.l) 


48 







For c = 1, all the entries of the Kac table are corner entries, so J>. s is completely reducible for all 
r, s € Z. In particular, we have 

Ai,3 Ai i5 Aij 

© • 0 • 0 • 0 ••• ( C - 2 ) 

(A\ not in Kac table). (C.3) 

As discussed in Section 5.5, the dimer Fock space (C 2 )®°° splits as the direct sum of subspaces 
E x . The corresponding bases are given by states of the form (5.86) with j — i = x for n odd and 
j — i = x + ^ for n even. Using (5.85) at k = 0 and the form of the basis, we find that the character 
over E x is a Verma module character, 


F\a 


Ex 


A 


1,1 


A a 


ch [E x ](q) 



nr=i(wr 


(C.4) 


where c = 1 and A 2 ^ is defined in (5.88). Crucially, for x € ^Z\ {0}, there are no integers r and s such 
that A 2 f. = Tj- equals A r . )S = ^ . In these cases, A 2 ^ is not in the c = 1 Kac table, and the character 

of E x is irreducible, thus proving the part of Theorem 2 pertaining to x 0. 

What remains is to prove the case x = 0, for which n is odd and the fermionic labels take half¬ 
integer values. We proceed by showing that E° — Eq by explicitly constructing the isomorphism. In 
fact, the proof presented below works in full generality, showing that E x — E\ for A = —x. As in 
Appendix B, we define two maps 


0 : H -> End((C 2 )®°°) and \) : E x -> E x , (A = -x), (C.5) 

defined as 

fl(l) = I) fl( a fc) = ~ Q-q-\-k^q ~ 5^,0 ^n,0 mod 2 

g^O g>0 

and 

• • • a Z\ x )) = 0(«u) m sK) n2 ■ • • 0K) n 1^)- 

We find that 

[g(a k ),g(a e )} = kd k+£fi , q(oo) = -X, g(a k )\£l x ) = 0 (k > 0), (C.8) 

implying that g is an H-homomorphism and E x a highest-weight H-module. This is almost the same 
as what was found in Appendix B. One key distinction is that the constant terms in fl(ao) an d 0 («o) 
are different. At the level of H, these constants are irrelevant, as any H-homomorphism g is defined up 
to an additive constant in g(ao). The choice of this constant is, however, important in the Virasoro 
picture. Here, we find that A^ F = A 2 ^ for A = —x, and 


(C. 6 ) 

(C.7) 


[Lf. ,a,q] — q a k + q , [L fc ,g(a g )] — qQ(a k +q)- 


(C.9) 


We note that modifying the constant in g(ao) would alter the second relation. We now have all the 
ingredients to show the next proposition. 


Proposition 3 The map b is bijective and it is an isomorphism of Virasoro modules. 

The proofs of the two statements in Proposition 3 are identical to those of Propositions 1 and 2 in 
Appendix B. It follows that E x ~ T— x and that the module structures of E x can be read off from (C.2) 
and (C.3), thus completing the proof of Theorem 2. 
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